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Preface 


This  research  topic  was  suggested  by  Dr.  Edward  J. 
Dunne,  Adjunct  Professor  in  the  Operational  Sciences 
Department,  School  of  Engineering,  Air  Force  Institute 
of  Technology,  Wright-Patterson  AFB,  OH. 

I  have  found  the  field  of  runway  closure  to  be  very 
relevant  today,  both  .in  the  probability  theory  and  in 
the  renewed. interest  in  attack  strategies.  Hopefully, 
the  development . of  a  computer  program  to  aid  in  planning 
runway  attacks  will  provide  further  insights  and 
experience  for  planners. 

I  am  grateful  to  my  thesis  advisor.  Dr.  Dunne,  for 
the  assistance,  encouragement  and  support  he  provided 
during  this  effort.  I  am  also  indebted  to  Lt.  Col. 

James  N.  Bexfield  for  his  understanding  and  technical 
advice. 
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Abstract 

A  computer  program  to  aid  the  planning  of  runway 
attacks  is  developed.  Conventions " ,  individually  targeted 
weapons  are  used  against  non-reinforced  concrete  runways. 
The  program  has  two  main  sections.  The  first  section 
evaluates  any  attack  strategy,  based  on  independent  cuts 
along  the  runwav,  with  each  cut  specified  in  terms  of 
number  of  aim  points,  number  of  weapons  per  aim  point,  and 
aim  point  locations.  The  second  section  searches  for  the 
"best"  strategv  which  uses  the  least  number  of  weapons  to 
achieve  an  overall  probability  of  runway  closure  equal  to 
or  greater  than  a  user  specified  level. 

The  program  operates  in  three  modes.  The  mode  1 
program  returns  the  fewest  number  of  weapons  and  the  "best" 
strategy  in  order  to  meet  or  exceed  a  user  defined  level  of 
runway  closure.  Mode  2  allows  the  user  to  specify  a  fixed 
number  of  weapons  instead  of  a  level  of  runway  closure,  and 
the  program  returns  the  highest  probability  of  runway 
closure  and  the  "best"  strategy  to  use  with  the  fixed 
number  of  weapons.  Finally,  mode  3  allows  the  user  to 
completely  specify  a  strategy  in  terms  of  number  of  cuts, 
cut  locations,  number  of  aim  points  per  cut,  number  of 
weapons  per  aim  point  and  locations  of  the  aim  point,  and 
the  program  returns  the  expected  probability  of  runway 
closure  for  the  user  defined  strategy. 
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I .  Introduction 


Background 

The  airfield  complex,  normally  attacked  in  hopes  of 
reducing  the  enemy1 s  air  capability,  is  a  popular  target 
for  military  strategists.  Of  primary  importance  are  the 
aircraft  because  of. thier  delicate  structure,  fuel  content, 
and  high-explosive  weapons.  Two  factors,  however,  make 
aircraft  less  opportune  targets  than  airfields.  First, 
aircraft  parked  in  shelters  are  difficult  to  destroy  since 
the  attacker  must  breach  or  substantially  damage  the  shelter 
in  order  to  reach  the  aircraft.  Special  weapons,  whose 
accuracy  in  some  cases  demand  special  aerial  maneuvers,  make 
attacking  aircraft  in  shelters  vulnerable  to  ground  defense 
systems.  Secondly,  increases  in  aircraft  performance  make 
aircraft  ^ess  opportune  targets.  Today,  attacking  aircraft 
can  strike,  return  to  base,  rearm,  and'  be  over  the  target 
again  in  an  extremely  short  time,  often  shorter  than  that 
required  to  repair  the  damaged  runway  before  defending 
interceptors  may  be  launched.  Thus,  opposition  fighters 
can  be  effectively  eliminated  from  the  first  stages  of  battle 
without  being  physically  damaged  by  denying  them  a  launch 
and  recovery  surface.  The  concept  that  runways  .  .  .  are 
poor  targets  is  no  longer  valid  in  the  light  of  the  pace  of 
modern  warfare  and  the  capabilities  of  aircraft.  (Ref  11) 

Currently,  development  of  ’’better”  weapons  for  runway 
destruction  is  the  emphasis  in  weapons  development.  An 
example  is  the  French-made  Durandal  runway  attack  weapon. 


which  is  a  parachute-retarded,  rocket-boosted,  concrete- 
penetration  bomb  designed  to  make  craters  of  runway  surface 
It  is  now  a  part  of  the  Foreign  Weapons  Evaluation  Program, 
and  if  further  testing  is  successful,  deliveries  will  be 
scheduled  for  mid-19 83.  (Ref  2) 


Problem  Statement 


The  following  scenario  forms  the  basis  for  this 
research.  How  should  a  wing  commander  allocate  his 
resources  to  effectively  close  a  runway?  The  runway  is  a 
rectangular  non-reinf orced  concrete  area  for  the  purpose 
of  launching  and  recovering  fixed-wing  aircraft.  The  run¬ 
way  is  closed  when  sufficient  space  to  launch  his  aircraft 
does  not  exist.  This  space  is  generally  thought  of  as  a 
rectangle  and  will  be  referred  to  as  the  minimum  launch 
window  (MLW). 

The  commander  has  available  only  one  type  of  indivi¬ 
dually  targeted  conventional  runway  munition  for  the 
mission.  The  attack  consists  of  cutting  the  runway  at 
specified  locations  along  the  length  of  the  runway.  Each 
cut  is  a  set  of  aim  points  located  along  a  line  perpen¬ 
dicular  to  the  runway  length  (see  Figure  1).  A  strategy 

\ 

is  a  set  of  completely  specified  cut  locations,  aim  points 
per  cut  and  number  of  weapons  per  aim  point. 


Figure  1  .  Attacking  Strategy 
For  reference,  the  origin  of  the  x-y  plane  will  be  at  the 


southwest  corner  of  the  runway  (see  Figure  1)  and  the 
length  of  the  runway  will  be  on  the  x-axis.  The  cut  loca¬ 
tions  are 'then  the  distances  along  the  x-axis.  The  aim 
point  locations  are  denoted  by  the  cut  location  and  the 
distance  from  the  edge  of  the  runway. 

A  review  of  available  literature,  discussed  in  Chapter 
II,  shows  that  the  targeting  strategies  for  runway  destruc¬ 
tion  has  not  kept  up  with  the  development  in  weapons. 

There  are  only  a  few  models  to  evaluate  targeting  strate¬ 
gies  and  none  which  will  find  a  "best”  strategy.  Thus, 
the  commander,  or  any  military  planner,  has  to  rely  on 
experience  or  insight  to  come  up  with  a  targeting  strategy 
before  using  one  of  the  existing  models. 

This  research  effort  develops  a  computer  program, 
described  in  Chapter  V,  which  will  find  a  "best”  targeting 
strategy  for  a  given  problem  without  requiring  the  commander 
to  come  up  with  a  targeting  strategy  himself.  This  computer 
program  will  aid  the  commander  by  answering  these  three 
questions : 

1.  How  many  weapons  does  it  take  to  obtain  a  certain 
level  of  runway  closure? 

2.  What  is  the  highest  level  of  runway  closure  that 
can  be  obtained  with  a  fixed  amount  of  weapons? 

3.  What  level  of  runway  closure  can  be  obtained  with 
an  independently  specified  strategy? 


Objectives 

The  objective  of  this  research. is  to  develop  a 
computer  program . which  will  answer  the  above  questions. 
Sub-objectives  are: 

1 .  Find  or  develop  a  math  model  to  evaluate  a 
strategy. 

a.  Model  must  be  capable  to  evaluate  any 

strategy. 

■  b.  Model  must  be  relatively  fast  in  terms  of 
computer  CPU  time. 

2.  Find  or  develop  a  search  algorithm  to  identify 


the  "best”  strategy  for  a  given  mission 


In  the  development  of  the  search  to  find  "best” 
strategies  or  evaluate  any  given  strategy,  the  following 
assumptions  are  made: 

1.  The  weapons  of  any  strategy  are  independent, 
individually  targeted. 

2.  The  error  distribution  of  any  weapon  is  circular 
normal  centered  at  the  aim  point  with  a  specified  CEP 
value. 

3.  Only  one  type  of  weapon  will  be  used  to  close  the 
runway  for  a  given  strategy. 

4«  The  survivability  of  all  weapons  is  one. 

5.  The  reliability  of  all  weapons  is  one. 

6.  There  is  no  bomb  damage  assessment  capability  to 
update  the  targeting  strategy. 

7.  The  minimum  launch  window  (minimum  length/ 
minimum  width  needed  for  takeoff)  must  be  composed  of 
actual  runway. 

8.  If  a  minimum  launch  window  exists,  it  is  assumed 
to  be  accessible  to  an  aircraft  for  takeoff. 

Scope  and  Limitations 

This  model  considers  only  conventional  weapons  in 
destroying  non-reinforced  concrete  runways. 

General  Procedure 

Basically  three  procedures  are  developed  in  this 
research. 


1.  Starting  from  the  independent  strategy,  the  prob- 


lem  is  to  calculate  the  probability  of  runway  closure  when 
the  strategy  is  completely  specified.  The  procedure 
developed  here  relies  on  either  a  discrete  approximation 
or  a  monte-carlo  simulation.  The  discrete  approximation 
is  faster  (on  the  computer)  for  problems  where  the  effective 
minimum  takeoff  width  is  approximately  half  or  more  of  the 
effective  runway  width.  If  not,  then  the  number  of  calcu¬ 
lations  necessary  for  the  discrete  approximation  makes  it 
less  efficient  than  the  monte-carlo  simulation,  which  is  a 
constant  run-time  routine. 

2.  The  next  problem  is  the  highest  probability  of 
runway  closure  for  a  given  number  of  weapons.  This 
problem  requires  a  search  procedure  that  will  consider  the 
possible  combinations  of  weapons  and  aim  points  and  aim 
point  locations  to  obtain  the  highest  probability  of  run¬ 
way  closure.  The  various  combinations  of  number  of 
weapons  per  aim  point  and  aim  point  locations  are  evaluated 
by  procedure  1  to  find  the  strategy  which  yields  the  high¬ 
est  probability  of  runway  closure. 

3.  The  last  problem  requires  the  level  of  runway 
closure  to  meet  or  exceed  a  predetermined  level  of  runway 
closure.  The  total  number  of  weapons  will  be  varied 
between  a  minimum  and  a  maximum  amount.  Each  number  of 
weapons  is  individually  searched  by  procedure  2  to  obtain 
the  highest  probability  and  stops  for  the  smallest  number 
of  weapons  that  meets  or  exceeds  the  predetermined  level 
of  runway  closure.  This  number  is  reported  as  the  ’’best" 


number  of  weapons  needed  to  obtain  the  predetermined  level. 


Organization 

The  first  phase  involves  finding  a  computer  routine 
that  can  evaluate  a  given  strategy  to  answer  the  following 
question:  What  level  of  runway  closure  can  be  obtained 
with  an  independently  specified  strategy?  Various  existing 
modeling  approaches  are  reviewed  in  Chapter  III.  In 
Chapter  III,  the  specific  computer  routine  will  be  described 
which  was  chosen,  based  on  both  speed  and  accuracy.  A 
relatively  fast  computer  routine  is  necessary  because  of 
its  repeated  use  in  the  search  algorithm. 

The  second  phase  involves  selecting  a  search  routine 
that  will  find  the  "best"  strategy  for  a  given  number  of 
weapons  to  answer  the  question:  What  is  the  highest 
level  of  runway  closure  that  can  be  obtained  with  a  fixed 
amount  of  weapons?  The  search  routine,  described  in 
Chapter  IV,  searches  among  the  combinations  of  number  of 
weapons  per  aim  point  and  aim  point  locations  to  find  the 
strategy  that  gives  the  highest  level  of  runway  closure. 

Finally,  the  third  phase,  involves  writing  the  overall 
program  that  will  find  the  "best"  strategy  to  meet  or 
exceed  a  predetermined  level  of  runway  closure  to  answer 
the  question:  How  many  weapons  does  it  take  to  obtain  a 
certain  level  of  runway  closure?  This  program,  described 
in  Chapter  V,  searches  for  the  minimum  amount  of  weapons 
necessary  to  achieve  a  certain  level  of  runway  closure. 


II.  Background 


This  chapter  is  a  review  of  the  currently  available 
literature  pertaining  to  probability  models.  Three  areas 
of  interest  are  discussed.  The  first  area  is  the  broad 
subject  of  coverage  problems  where  the  probability  of 
destroying  a  fixed  target  is  given.  Next,  is  a  review  of 
models  that  evaluate  a  particular  runway  closure  strategy. 
The  third  is  a  review  of  the  one  model  found  that  incor¬ 
porates  a  searching  scheme  with  an  evaluation  method, 
resulting  in  a  targeting  strategy. 

Coverage  Problems 

The  evaluation  of  targeting  strategy  can  be  logically 
developed  by  first  considering  the  impact  of  a  single 
weapon  on  a  fixed  rectangular  target  (e.g.,  runway). 

This  type  of  research  comes  under  the  general  heading  of 

coverage  problems.  One  very  helpful  source  in  this  area 

\ 

is  the  "Survey  of  Coverage  Problems  Associated  with  Point 
'and  Area  Targets"  by  A-  R.  Eckler  (Ref  3).  Jaiswan  and 
Sengal  looked  at  the  problem  of  the  expected  damage  area 
for  stick  and  triangular  pattern  bombing  (Ref  7).  All 
these  results  consider  only  fixed  targets;  however,  the 
problem  of  runway  closure  must  consider  the  minimum 
launch  window  (i.e.,  the  least  amount  of  runway  surface 
necessary  for  aircraft  launch,  as  the  actual  "target".  In 
other  words,  the  objective  of  bomb  damage  is  to  deny  an 
area  equal  to  the  minimum  launch  window  anywhere  on  the 
runway.  Thus,  just  computing  the  expected  area  damaged 


without  regard  to  where  the  damage  occurs  is  not  sufficient 
for  calculating  the  probability  of  runway  closure. 

Probability  of  Closure  Models 

The  models  that  can  explicitly  compute  the  probability 
of  runway  closure  are:  AIDA,  TSARINA,  AHAB,  RUNW  and  Manz. 
AIDA,  TSARINA,  AHAB  and  RUNW  are  all  monte-carlo  simulations. 

AIDA  (Airbase  Damage  Assessment ). program  (Ref  4)> 
which  is  the  principle  model  used  by  USAF/Studies  and 
Analysis,  is  a  large  scale  FORTRAN  IV  model  that  considers 
the  whole  airfield  complex  with  up  to  250  separate  targets 
in  the  complex.  Targeting  flexibility  includes  250  sepa¬ 
rate  targets,  such  as  parked  aircraft,  POL  centers,  hangars, 
etc.  Each  target  can  be  assigned  to  a  separate  vulnerabil- 
ity  class;  hard  targets  such  as  hardened  hangars  may  be 
one  vulnerability  class  while  aircraft  in  the  open  are  in 
another.  Up  to  20  different  vulnerability  classes  may  be 
identified. 

Each  weapons  delivery  pass  is  described  by  weapon 
type  and  weapon  releaae  parameters  (such  as  heading, 
altitude,  air  speed,  dive  angle,  and  stick  length).  Up 
to  10  different  weapon  types  may  be  used  in  combination 
with  up  to  50  separate  weapon  delivery  passes.  Also 
considered  are  the  probability  of  arriving  at  the  target, 
attrition  due  to  enemy  fire,  aiming  accuracy  in  terms  of 
REP  and  DEP  and  ballistic  dispersion. 

With  regards  to  runway  targets,  only  point-impact 
'■■v  weapons  are  considered.  The  effective  miss  distance  is 


10 


equal  to  the  damage  radius.  Up  to  250  hits  may  be  stored 
for  each  runway.  Should  the  runway  be  closed,  it  identifies 
the  minimum  number  of  craters  that  have  to  be  repaired 
before  the  runway,  can  be  reopened.  An  approximate  computer 
plot  of  impact  points  is  available  upon  request. 

The  probability  of  runway  closure  is  calculated  by  the 
proportion  of  times  the  runway  was  closed  during  the  simu¬ 
lation.  The  runway  is  considered  closed  if  no  rectangle 
equal  to  the  minimum  launch  window  exists  anywhere  on  the 
runway  without  a  crater.  AIDA  determines  this  by  starting 
at  one  corner  of  the  runway,  positioning  a  rectangle  equal 
to  the  minimum  launch  window  on.  it  and  then  observing 
whether  or  not  the  rectangle  contains  a  crater.  If  it 
does  contain  a  crater,  the  rectangle  is  then  shifted  five 
feet  along  the  width  of  the  runway  and  is  again  checked 
for  craters.  If  an  open  rectangle  does  not  exist  (each 
rectangle  contains-  one  or  more  craters)  when  the  procedure 
reaches  the  other  end  of  the  runway  width,  the  rectangle 
is  then  shifted  250  feet  along  the  length  of  the  runway, 
positioned  at  one  edge  of  the  runway  width,  and  the  search 
continues  until  either  an  open  rectangle  is  found  (the 
entire  runway  is  open)  or  the  procedure  reaches  the  oppo¬ 
site  corner  of  the  runway  and  no  open  rectangle  is  found 
(the  entire  runway  is  closed). 

TSARINA  (Ref  5)  is  a  version  of  AIDA  that  makes  the 
AIDA  model  compatible  with  the  TSAR  (Theater  Simulation  of 
Airbase  Resources)  model,  thereby  building  a  very  large 


scale  model 


AHAB  (Attacking  Hardened  Air  Bases)  by  RAND  Corporation 
(Ref  10)  is  an  interactive  computer  simulation  designed  to 
aid  decision  makers  in  finding  an  airbase  attack  plan  which 
maximizes  their  value  function.  The  value  function  must  do 
three  things: 

1 .  Reflect  the  natural  ordering  of  preferences  -  in 
fact,  rank  losing  two  of  your  aircraft  as  superior  to 
losing  three  of  them. 

2.  Express  various  components  of  an  outcome  in 
comparable  terms.  Is  losing  two  of  your  own  planes  to 
destroy  four  of  the  enemy's  better  or  worse  than  losing 
five  to  destroy  eight? 

3.  Distinguish  among  uncertain  outcomes.  For 

example,  although  we  may  know  that  there  is  a  fifty-fifty 
chance  that  an  attacking  aircraft  will  destroy  a  hangarette, 
we  cannot  know  how,  many  hangarettes  will  be  destroyed  ( 

during  a  particular  mission,  nor  can  we  be  certain  of  how 
many  aircraft  will  be  in  those  hangarettes.  The  value 
function  places  values  on  the  probability  distributions  of 
outcomes.  We  should  be  able  to  decide  whether  a  50  percent 
chance  of  losing  six  aircraft  and  a  50  percent  chance  of 
losing  none,  is  better  than,  worse  than,  or  the  same  as 

100  percent  chance  of  losing  three,  all  other  things  being 
equal. 

The  von  Neuman-Morgenstern  utility  function  is  the 
value  function  used  in  this  model.  This  function  reflects 
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not  only  the  appropriate  values  or  utilities  of  outcomes, 
but  a  characterization  of  the  decision  maker’s  attitudes 
toward  risk. 

The  value  function  is  linear  and  has  four  arguments: 
the  number  of  enemy  aircraft  destroyed,  the  number  of 
aircraft  shelters  or  hangarettes  destroyed,  the  number  of 
hours  the  runway  complex  is  closed,  and  the  number  of 
friendly  aircraft  lost  in  the  attack. 

AHAB  is  an  interactive  monte-carlo  simulation  written 
in  JOSS  (a  RAND  language).  The  target  complex  consists  of 
open  aircraft,  hangarettes  and  runways.  Only  one  weapon 
type  is  considered  for  a  problem,  and  the  attack  is  assumed 
to  come  across  the  runway  perpendicular  to  the  runway  at 
evenly  spaced  ’’cuts"  along  the  runway. 

The  calculation  of  the  probability  of  runway  closure 
is  to  look  for  a  gap  of  50+d  feet,  assuming  a  minimum  take¬ 
off  width  of  50  feet,  along  any  of  tiie  ’'cuts"  where  d  is 
equal  to  two  times  the  damage  radius  of  the  weapon.  If  a 
gap  of  50+d  is  found  along  any  "cut"  then  the  runway  is 
open,  if  not,  then  the  runway  is  closed. 

RUNW  (Ref  12)  is  a  NATO  program  specifically  designed 
for  the  calculation  of  runway  interdiction  probability. 

RUNW  is  a  monte-carlo  simulation  model  that  computes  only 
the  probability  of  closure  for  one  runway.  The  target  is 
a  rectangular  runway  and  the  minimum  launch  window  is  a 

rectangle  (minimum  launch  length  by  minimum  launch  width). 
An  attack  is  described  by  the  weapon  type,  attack  heading 


(from  parallel  to  the  runway  to  perpendicular  to  the  runway), 
whether  the. weapons  are  released  individually,  in. salvos  or 
in  sticks,  the  aiming  error  (normally  distributed),  the 
weapon  dispersion  error  (normally  or  uniformly  distributed). 

The  probability  of  runway  closure  for  RUNW  is  the  pro¬ 
portion  of  times  the  runway  is  closed  over  the  number  of 
iterations  used.  The  procedure  for  finding  an  open  minimum 
launch  window  first  orders  the  impact  points  on  the  runway 
along  the  x-axis  (length)  then  looks  at  the  subsequent 
spaces  between  these  ordered  points.  If  any  space  is  greater 
than  the  minimum  launch  length,  then  the  runway  is  open.  If 
there  are  no  spaces  greater  than  the  minimum  launch  length, 
it  then  looks  at  spaces  with  one  point  in  between.  If  there 
is  a  space  greater  than  the  minimum  launch  length  now,  the 
lateral  spacing  between  the  runway  edge  and  the  bomb  is 
checked  against  the  minimum  launch  width,  and  if  greater, 
then  the  runway  is  open.  The  search  continues  until  the 
case  when  four  interior  points  exist  in  the  spacing.  If  no 
minimum  launch  window  is  found,  then  the  runway  is  considered 
closed. 

The  model  by  Dr.  Manz  (Ref  9)  calculates  the  probabil¬ 
ity  of  runway  closure  for  cluster  munitions  that  are 
uniformly  distributed.  This  is  an  analytical  model  based 
on  the  binomial  distribution  of  submunitions  falling  in  a 
given  subarea  of  the  runway  and  a  differential  equation 
argument  to  compute  the  probability  of  runway  closure. 

The  probability  of  runway  closure  P(cJnq),  given  that 
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'-V‘. 


N  submunitions  are  aimed  at  the  runway,  is 
o 


a 

p(c|N0)  =  p(c|n)  p(n|n0) 


where  P(n|Nq)  is  the  probability  that  N  submunitions  will 
impact  on  the  runway  given  that  Nq  submunitions  were  aimed 
at  the  runway.  This  is  a  binomial  distribution  of  the 
submunitions  with  a  probability  Pq  for  any  ’one  submunition 
impacting  the  runway. 

P(c|n)  is  theprobability  of  runway  closure  given  that 
N  submunitions  have  impacted  the  runway.  This  probability 
can  be  evaluated  for  a  given  rectangle  on  the  runway.  Then 
the  rectangle  is  moved  a  distance,  dx,  in  the  x  direction, 
and  dy  in  the  y  direction.  The  probability  that  a  rectangle 
that  was  closed. will  be  opened  up  by  the  move  is  calculated 
by  solving  a  differential  equation. 

Runway  Closure  Strategy  Model  " 

Only  one  model  has  attempted  to  search  for  the  "best” 
strategy  in  terras  of  aim  point  locations  and  number  of 
weapons  per  aim  point.  The  "best”  strategy  is  defined  as 
the  strategy  which  will  obtain  the  required  level  of  runway 
closure  with  the  least. number  of  weapons.  This  model  is 
by  Gapt.  Pemberton  (Ref  11)  and  contains  two  parts. 

The  first  calculates  the  probability  of  runway  closure 
for  a  given  strategy  by  either  a  discrete  approximation  or 
a  monte-carlo  simulation.  The  technique  used  depends  on 
the  speed  of  the  calculation.  The  monte-carlo  simulation 
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requires  a  lot  of  computer  time  to  evaluate  any  strategy. 

The  discrete,  approximation  evaluation  time  depends  on  the 
dimensions  of  the  runway  and  the  MLW.  When  the  MLW's 
width  is  from  1/3  to  approximately  equal  to  the  runway's 
width,  the  discrete  approximation  requires  fewer  calcu¬ 
lations  and  is  faster  than  the  monte-carlo  simulation. 

When  the  MLW' s  width  decreases  in  proportion  to  the  runway's 
width,  the  number  of  calculations  increase,  and  when  the 
MLW's  width  is  approximately  equal  to  1/3  of  the  runway's 
width,  the  number  of  calculations  required  in  the  discrete 
approximation  make  the  calculating  time  greater  than  that 
required  by  the  monte-carlo  simulation.  Thus,  for  situ¬ 
ations  where  the  MLW's  width  is  less  than  1/3  of  the  run¬ 
way's  width,  the  program  uses  the  monte-carlo  simulation 
because  it  is  faster. 

The  second  part  is  a  search  routine  which  finds  the 
"best"  strategy.  The  routine  systematically  searches  over 
the  combinations  of  aim  point  locations  and  number  of 
weapons  per  aim  point  by  first  selecting  the  "best"  aim 
point  locations  then  allocating  weapons  to  these  aim  points 
until  the  required  level  of  runway  closure  is  obtained. 
Comparison  of  Models 

From  this  research  it  was  concluded  that  the  model  by 
Pemberton  was  the  closest  to  meeting  the  requirements  of 
the  problem.  However,  the  Pemberton  model  was  considered 
to  have  weaknesses  in  the  discrete  approximation  section 
and  needed  a  better  search  routine.  The  objectives  of 
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this  effort  are  to  improve  the  Pemberton  model  and  apply 
the  model  in  a  program  which  could  be  a  useful  decision 
aid  to  a  military  planner. 


Ill .  Evaluation  of  a  Given  Strategy 


Expansion  of  the  problen  to  throwing  a  point  at  an  extended 
target 

When  considering  weapons  with  a  damage  radius  greater 
than  zero,  the  impact  point  of  the  weapon  need  not  coin¬ 
cide  with  the  target  to  damage  it.  While  the  weapon  may 
miss  the  target  by  an  amount  equal  to  the  radius,  it  can 
still  effectively  damage  the  target.  On  the  other  hand, 
by  expanding  the  target’s  boundaries  by  an  amount  equal 
to  the  damage  radius  all  around  and  considering  the  weapon 
as  a  point  weapon,  damage  to  the  target  could  then  be 
determined  by  observing  whether  or  not  the  impact  point  is 
within  the  extended. boundaries  (see  Figure  2). 


r* 

-impact  point 

— - 

r 

Target 

Target 

J 

Extended  Target 

Figure  2.  Extended  Target  Boundary 
In  this  research,  runway  and  minimum  launch  window 
dimensions  are  converted  into  extended  dimensions  for  the 
probability  calculations.  For  example,  the  extended 
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minimum  takeoff  width  (ETOWID)  is  equal  to  the  minimum 
takeoff  width  (TOWID)  plus  twice  the  damage  radius  (DAMRAD). 
Reduction  of  the  two  dimensional  problem  to  a  one 
dimensional  problem 

The  ideal  approach. to  the  problem  of  computing 
runway  closure  probability  involves  computing  the  probability 
that  an  open  space,  in  both  length,  and  width  does  not  exist 
anywhere  within  the  runway  boundaries.  For  an  overall 
approach,  one  can. order  all  of  the  weapons  that  impact  on 
the  runway  in  increasing  order,  from  the  origin,  see  figure  3. 


Figure  3.  Ideal  approach  for  probability  of  closure 


Under  this  arrangement,  rectangular,  spaces  are  formed  by 
either  four  points  or  by  the  runway  boundaries  and  impact 
points.  These  rectangles  are  then  compared  to  the  minimum 
launch  window.  If  any  rectangle  is  larger  than  the  MLW, 
then  the  runway  is  open.  The  probability  that  any 
rectangle  is.  larger  than  the  MLW  depends  on  the  joint 
probability  distribution  of  all  impact  points. 

This  two-dimensional  problem  can  be  reduced  to  a  one 
dimensional  problem,  since  runways  are  many  times  longer 


than  they  are  wide.  The  closing  of  a  given  MLW  by  denying 

it3  length  is  dependent  on  closing  of  that  MLW  by  denying 

its  width.  See  figure  4.  Runway  closure  can  be 

accomplished  if  we  have  one  "cut”  per  MLW,  where  a  cut  is 

defined  to  be  a  strategy  that  attempts  to  deny  a  space 

along  a  line  perpendicular  to  the  runway  length  equal  to 

the  MLW  width.  Then  the  probability  that  a  MLW  is  closed  is 

* 

equal  to  the  probability  of  cut.  Conversely,  the  probability 
that  a  MLW  is  open  is  equal  to  the  probability  of  not  cutting 
the  runway  at  that  location.  This  assumes  that,  given  a 
successful  cut,  no  weapon  extends  beyond  the  length  of  the 
MLW.  This  can  be  virtually  assured  by  aiming  the  cut  at 
least  three  sigma  (see  appendix  B  for  the  relationship 
between  CEP  and  sigma)  away  from  either  edge  of  the  MLW. 

Thus,  the  probability  of  closure  is  then  the  product 

of  the  probabilities  of  cut  of  the  individual  cuts. 

\ 

Arrangement  of  cuts  along  the  runway  and  the  minimum 
launch  lengths  are  shown  in  figure  4. 


3CT  1  3<T 

y5  '  3C 

cut  1 

cut  2 

Figure  4.  Individual  Cut  Locations 


The  runway  is  open  if  any  cut  fails  to  deny  the  minimum 
launch  width.  For  example,  if  cut  1  fails  then  the  length 
from  the  left  edge  of  the  runway  to  cut  2  is  open:  Since 
there  is  a  gap  greater  than  the  minimum  launch  width  in 
cut  1  there  is  at  least  one  open  minimum  launch  window 
(MLW1).  If  cut  2  fails,  then  a  minimum  launch  window  can 
be  positioned  in  the  space  in  cut  2  with  sufficient  length 
between  cuts  1  and  3;  thus  the  runway  is  open. 

The  runway  is  closed  only  if  all  cuts  succeed.  The 
number  of  cuts  required  is  calculated  by 

NCUTS  =  '( ERUNLE  -  6(STADEV))1  ^ 

.  SRL  -  6(STADEV )  J 

where  NCUTS,  the  number  of  cuts  required,  is  the  greatest 
integer  less  than  the  quantity  in  the  brackets.  STADEV  is 
the  standard  deviation  of  the  impact  point  distribution 
(sigma)  derived  in  appendix  B.  SRL  is  the  shortest  runway 
length  for  takeoff  derived  in  appendix  C.  Thus,  the 
probability  of  runway  closure  is 


NCUTS 

’cl  =  7T  poi 

CL  i=1  1 


where  PC^  is  the  probability  of  cut  for  location  i.  This 
method  assumes  that  sigma  is  less  than  ETOLEN/6,  where 
ETOLEN  is  the  effective  takeoff  length,  because  of  the 
requirement  for  independence  of  the  cuts  in  this  case. 

We  can  now.  concentrate  on  the  probability  of  cut  (PC). 
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This  is  the  probability  of  denying  (along  the  width)  a 
space  equal  to  the  minimum  launch  width.  Three  methods  of 
calculating  the  probability  of  cut  were  examined:  an  order 
statistics  approach,  a  discrete  approximation  approach  which 
approximates  the  continuous  nature  of  successive  minimum 
widths  by  discrete  steps,  and  monte-carlo  simulation. 

Order  Statistics  Approach 

The  problem  is  to  calculate  the  probability  that  the 
largest  space  along  a  line  perpendicular  to  the  length  of 
the  runway  is  less  than  a  certain  minimum  launch  window's 
width.  the  largest  space  is  the  max  |x^  +  ^j  -  x(-jjj 

where  x(j_)  are  the  order  statistics  of  the  impact  points 
from  one  edge  of  the  runway,  with  x^gj  =  0  and  x^n+1 ^  = 
runway  width. 

If  we  assume  a  uniform  distribution  across  the  runway 
and  assume  that  the  runway  width  is  1 ,  then  the  probability 
that  £  is  greater  than  v  (the  minimum  launch  window's  width) 
has  been  derived  by  David  (Ref  1),  and  is  given  by: 

Pr(£>vj  =  n(1-v)n“1  -Q)  (1-2v)n_1  +  ••• 

(4) 


n  =  number  of  spaces  =  number  of  weapons  +  1 
where  the  series,  continues  as  long  as  (1-iv)>0.  (Ref  1:81) 
The  problem  in  this  research  is  to  apply  this  procedure 
not  to  independent  identically  distributed  uniform 
distributions,  but  to  independent  non-identical  normally 
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distributed  bombing  distributions.  Unfortunately,  when 
considering  the  normal  distribution,  the  order  statistics 
approach  gets  more  complex.  When  the  order  statistics  of 
the  joint  normal  distribution  are  superimposed  on  a  fixed 
range  (runway  width)  this  introduces  two  more  order 
statistics  not  from  the  previous  joint  distribution.  The 
next  step  is  to  take  the  difference  between  the  order 
statistics  begining  from  one  end  of  the  runway  width  to 
the  other  end..  These  differences  are  then  ordered  and  we 
are  interested  in  the  largest  ordered  space  statistic.  The 
distribution  of  this  order  statistic  could  not  be  found, 
and  derivation,  though  perhaps  possible,  was  beyond  the 
scope  of  this  research  (Ref  6). 

Since  an  exact  solution  of  the  problem  could  not  be 
found,  the  problem  was  divided  into  discrete  subproblems 
discussed  in  the  next  section. 

Discrete  Approximation 

This  approach  yields  an  analytic  solution  based  on  the 
set  theory  concept  of  unions  and  intersections  of  non- 
independent  subevents  to  express  a  certain  event.  The 
event  of  interest  is  the  runway  cut  denoted  (RC).  The 
runway  is  not  cut  (i.e.,  open)  (denoted  RC)  if  there  exists 
a  space  equal  to  or  greater  than  the  minimum  launch 
window's  width.  The  runway  is  thought  of  as  containing 
some  number  of  discrete  overlapping  minimum  launch  window 


widths,  illustrated  in  figure  5.  These  overlapping 
minimum  launch  widths  are  called  sections  (s)  and  there 


are  NSPW  (number  of  sections  per  runway  width)  sections. 


NSPW  = 


The  s^'s  are  the  discrete  minimum  launch  window 


widths  that  are  to  be  evaluated.  The  event- 


closed 


occurs  when  a  weapon  lands  within  section  s^.  The 
event js^  openj occurs  when  all  weapons  miss  section  s^. 
NSPW  is  the  number  of  sections  s.^  per  runway  width. 

Figure  5.  Relation  between  s^’s  and  runway  width 
The  event  RC  is  equivalent  to  the  event^all  s^  are  closedj, 
where  the  event  £s^  closedj is  when  the  section  sustains 
some  damage  (at  least  one  impact  point  is  within  s^)  and 
the  event  £s^  openj is  when  the  section  does  not  sustain 
any  damage. 


NSPW 


RC  = 


a  n  / 

=  n  Is. 

i=i  1 


closed 


The  events  Js^  openj  are  not  independent  so  the  probability 
of  RC  denoted  (PC)  cannot  be  expressed  as  the  product  of  the 
individual  probabilities  of  s^.  The  event  RC  is  equivalent 
to  the  event ^at  least  one  s^  is  openj. 
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nu  =  \j  |s^  open^ 

So,  the  probability  that  the  runway  is  not  cut  PC, 

(Pr  {rc})  can  be  expressed  by  the  probability  multiplication 
law  as: 


, _ •)  _  NSPW 

|RCj  =  PC  =  £  Pr  |s^  open] 


NSPW-1  NSPW 


i=1  j  =2 
i<j  , 


^SPW 

,4,  Pr  W,sj 


open 


NSPW- 2  NSPW-1  NSPW  , 

+  ±5,  ^2  k?3  Prtwsk  openi 


i<j<k 

,  1 \ NSPW-1  p  J  o 
(  ~1 )  Pr  ]^s^ ,  S2» 


’  SNSPW  °Pen} 
(Ref  8:33) 


NSPW 

There  are  2  -1  terms  in  this  expression.  In  Pemberton's 

UNION  routine  (Ref  11:20)  all  2^^-1  terms  are  calculated. 
This  routine  can  be  made  more  efficient  by  observing  that 
many  of  these  terms  cancel  with  each  other.  The  event 
|s^,Sj  open]  occurs  when  no  weapon  damages  either  the  itJl 

X  T_ 

or  the  j  section.  If  s.  and  s.  overlap  (s.Os.^0)  then 

J  J 

the  event  [s^,Sj  open]  occurs  when  all  weapons  land  above 
the  upper  section  and/or  below  the  lower  section.  If, 
however,  there  exists  a  section  s^»  such  that  s^  is  wholly 
contained  in  fs.Os.},  then  the  event  { s. ,s.  open]  is 


equivalent  to  the  event  |si»sj»sjc  °penj.  In  like  manner 
all  events  that  include  s^,  s  j  and  wholly  contained  sections 
and  combinations  of  sections  in  js^ns^j  are  equivalent 
events  and  their  respective  probabilities  are  equal. 


It  can  be  shown  that  if  s^  and  overlap 


(s1flsNSP¥'!0)’  then  for  every  s^,  s^  (i/j^i+1)  there  exists 
equal  probabilities  with  opposite  signs  that  cancel  with 


each  other.  For  example,  if  NSPW  =  4*  and  (s.jfls^0)  then: 


PC 


=  it,  Pr(si 


open^  - 


>  )  PrSs.,s.  openr 

i=1  j  =2  1  J  J 


i<j 


2 

+  L 


i=1  jt=2  k“ 
i<  j<k 


Prfs.  ,  s.. ,  sv  open][ 


ri*  r  k 


-  Pr  js^s^s^.s^  open  j 


(8) 


=  Pr|s1  open] +  Pr{s2  °Pen]  +  Pr  js^  open]  +  Pr{s^  open] 

-  Pr|s^,S2  open]-  Pr[s^,s^  open]-  Pr^s^,s^  open] 

-  Prfs^s^  open]  -  Pr]s2»s^  open]-  Pr(s^»s^  open] 

+  Prfs^s^s^  open]+  Pr{s^,s2,s^  open] 

’  S4 

*  ^3  * 

Since  (s^Os^^0)  and  sections  s2and  s^  are  wholly  contained 
in  Is.fts.i  then  the  events: 


+  Pr  ,  s^ 
-  Pp{»1.b2 


js^s^  openj  =  open}  ) 

{S2,S4  °Pen]  —  {s2’s3,s4  °Penl 

are  equivalent  events  and  their  probabilities  are  equal. 

Each  probability  term  has  an  equal  probability  term  with  the 
opposite  sign  so  they  cancel  with  each  other  leaving  only 
these  terms: 

» 

Pr|s.|  open] +  Prjs2  open]  +  Pr  ^  open]  +  Pr{s^  open] 

-  Pr]s1,s2  open]-  Pr[s2,s^  open]-  Pr{sj,s^  open] 


_  NSPW  NSPW-1 

PC  =  ^  Pr{s.  open]-  ^  Pr^s.,3. 


open 


(10) 


The  probability  Pr^s.^  open]  is  equal  to  the  joint 
probability  that  every  weapon  impacts  outside  the  section  s^ 
For  example,  consider  two  weapons  and  xan  arbitrary  section 
s.  illustrated  in  figure  6. 


Figure  6.  Prjs^  open] 

The  probability  that  s^  is  open  from  weapon  1  is  the 
probability  that  weapon  1  impacts  above  3^  or  impacts 
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below  s^.  The  probability  that  s^  is  open  from  weapon  2  is 
the  probability  that  weapon  2  impacts  above  or  below  s^ . 

Since  weapons  1  and  2  are  independent,  that  is  weapon  Vs 
impact  has  no  influence  on  weapon  2's  impact,  the  probability 
that  s^  is  open  is  the  product  of  these  two  probabilities  or 

Pr{s.  open!  =  Pr  {weapon  1  impacts  above  or  below} 

L  1  J  \  J(11) 

•Pr  }  weapon  2  impacts  above  or  below} 

n 

In  general  Pr{s.  open}  =  7f  Pr {weapon  i  impacts  above  or 

L  J-  J  i=1 

below  section  i} ,  where  n  represents  the  total  number  of 
weapons  aimed  at  that  cut.  Because  s^  and  overlap,  the 

probability  Pr{s^,s^+^  open]  has  the  same  form  as  the 
probability  Pr^  open  j  but  the  impacts  are  now  above  and 
below  both  s^and  si+^ .  This  is  equivalent  to  the 
probability  that  each  weapon  impacts  above  the  upper  section 
and  impacts  below  the  lower  section.  In  general: 


Pr^s^s^-j  openj  =  7T  Prjweapon  i  impacts  above 


si+1  or  ^elow 


The  case  where  i=2  is  illustrated  in  figure  7. 


(12) 


Figure  7.  Pr{si»si+i 


NSPW  =  6 
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Figure  8.  sinsNSPW=^ 

If  and  s^Sp^  do  not  overlap,  then  some  of  the  pairs 
of  equivalent  events  identified  when  s^  and  sjjgpy  overlap 
do  not  exist.  These  pairs  have  been  identified  for  two 
cases. 

Case  I.  s-iflsNSPvf^’  s1^sNSPW-1^*  see  fiSure  8*  From  the 
figure  we  see  that  the  event  [si»smspw  °Pen]  is  n0't 
equivalent  to  the  event  £ s^ ,  S£»  s^gPW  °Pen|*  Thus  we  need 
two  extra  terms,  in  addition  to  the  terms  identified  in 
equation  10. 


NSPW  NSPW-1  . 

=  Z1  Pr(si  °Pen)  “  L1  Pr(si»si+1  °Pen5 

-  Pr{s1,sNSPW  °penj  +  Pr{s1’s2’sNSPW  °penj 


(13) 


where 


Pr{s1lSNSPW  °pen}  =  t  Prjweapon  i  impacts  above 
SNSPW*  between  sjjgpy  anc^  si  or  below 
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NSPW  =  6 


Figure  9.  s1OsNgpw_1=0 

Case  II.  si^sNSPW=^’  sinsNSP¥-1 =^’  SlnSNSP¥-2^’  see 
figure  9.  From  this  figure  we  see  that  the  event 

(s1’sNSP¥  °Pen^s  n0^  equivalent  to  the  event  fs-)  *  s2' SUSP¥ 
open].  The  event  (s-|»snsP¥-1  °Penj  no^  equivalent  to  the 

event  ( si  » s2»  sNSP¥-1  °Peni  and  the  event  {s2’sNSP¥  °Penj 
is  not  equivalent  to  the  event  j  Sg, s^» s^gpw  open],  and  the 

event  ^ s^ , sNgpw_1 , sNgp^  open j  is  not  equivalent  to  the 

event  (si » s2» sNSP¥-1 ’ sNSP¥  °Penj *  Thus  these  eight  terms 
need  to  be  added  to  the  terms  identified  in  equation  101 


*  sNSP¥  °Pen$  -  Pl1s2’SNSP¥  0pen^  "  Pris1’SNSP¥-1 

,S2’SNSP¥  openi  +  Prfs2,s3'sNSP¥  open3 

2*  NSP¥-1  p  3  +  Prf  s-|  »  s2>  sNSP¥-1  ’  SNSP¥  °Penj 
,S2,SNSP¥-1 ,SNSP¥  openJ 


Further  cases  can  be  examined  but  the  terms  needed 
in  1 

increase  as  2  -2  where  m  stands  for  the  largest  ra  such 


that  si^sNSPW-m^  true*  f°r  example,  if  the  largest 
non-overlapping  sections  are  s^  and  s^gp^  ^  then  m=3,  and 
the  number  of  extra  terms. to  be  added  to  the  terms  in 
equation  10  is: 

m  *  I  «* 

V  2  -2  =  2  +  6  +  14  =  22  extra  terms, 

i  =  1 

If  the  STEP  size,  which  is  the  distance  between  the 
lower  ends  of  two  adjacent  sections  is  equal  to  five  feet 
then  with  the  following  ranges  of  input  variables: 

RUNWID  -  (100,  150) 

TOWID  -  (  50,  150) 

W  -  (250,  2000) 

The  worst  case  is  where  sections  s^  and  s^gp^  ^  do  not  over¬ 
lap  which  are  the  conditions  for  case  II  discussed  above. 
Thus  using  this  modified  discrete  approximation  routine  the 
probability  of  cut  can  be  calculated  with  a  maximum  of 
2(NSPW)-1 +2+6+14  =  2(NSPW)-1 +22  =  2(NSPW)+21  terms.  The 

number  of  terms. necessary  for  ^emberton's  UNION  routine 

NSPW 

requires  calculating  all  2  -1  terms.  This  represents 

a  minimum  savings  of  2NSPW-1-2(NSPW)-21  =  2NSPW-2(NSPW)-22 

terms  that  need  not  be  calculated.  The  maximum  savings  is 

2NSPW-1-2(NSPW)+1  or  2NSPW-2(NSPW) .  For  example,  if 

NSPW  =  12,  the  new  approximation  routine  calculates  2(12) +21 

1  2 

or  -45  terms  and  Pemberton’s  routine  calculated  2  -1  or 

4095  terms.  The  savings  for  this  one  calculation  of  the 
probability  of  cut  is  4095-45  or  4050  terms  that  need  not 


be  calculated 


Problems  with  input  variables  that  exceed  the  above 
ranges  and  have  non-overlapping  sections  s^  and  sjjgpw_2  are 
approximated  with  a  monte-carlo  simulation. 


Figure  10.  PC  for  one  weapon 


A  special  case  of  the  discrete  approximation  is  when 
considering  the  probability  of  cut  when  only  one  weapon  is 
used.  If  s^  and  s^gpy  do  not  overlap,  then  the  probability 
of  cut  is  zero,  since  the  one  weapon  can  only  damage  one  of 
the  sections  leaving  the  other  section  undamaged.  The 
probability  of  cut  when  s^  and  do  overlap  is  the 

probability  that  the  weapon  impacts  in  the  intersection  of 
s^  and  Sjjspw  This  probability  is  just  the  area  under  the 
normal  error  distribution  of  the  weapon  for  the  intersection 
of  sections  s^  and  s^gp^  as  shown  in  figure  10.  The  weapon 
must  impact  in  the  shaded  area  in  order  to  cut  the  runway, 
thus  the  shaded  area  under  the  curve  is  the  probability  that 
the  weapon  will  cut  the  runway. 


Monte-Carlo  simulation 

The  monte-carlo  simulation  is  a  straight  forward 
approximation  of  the  real  world  probability  by  "attacking" 
the  runway  many  times  and  taking  the  proportion  of  times  the 
runway  was  cut.  Again,  runway  closure  is  the  product  of  the 


independent  probabilities  of  cut: 


Figure  11.  Simulated  Attack 

The  runway  is  represented  by  the  extended  runway  width 
(ERUNWI)  and  the  MLW  is  represented  by  the  extended  takeoff 
width  (ETOWID).  Runway  cut  is  simulated  by  generating 
weapon  impact  locations  (hits)  from  a  normal  random  number 
generator.  This  number  is  then  translated  by  multiplying 
it  by  the  standard  deviation  and  then  adding  the  mean  value 
for  that  weapon  distribution  (STADEV  and  aim  point  for  that 
weapon)  to  obtain  the  hit  locations.  After  all  weapons 
have  hit,  see  figure  11,  weapons  that  hit  below  or  above  the 
effective  runway  are  ignored.  If  none  of  the  weapons  hit 


the  effective  runway,  then  the  runway  is  not  cat.  If  there 
is  at  least  one  weapon  hit  on  the  effective  runway,  the 
distance  between  the  lower  edge  of  the  effective  runway  and 
the  nearest  hit  location  is  compared  to  the  effective 
takeoff  width,  if  the  distance  is  greater  then  the  runway 
is  not  cut.  Then  the  distance  between  the  upper  edge  of 
the  effective  runway  and  the  nearest  hit  location  is 
compared  to  the  effective  takeoff  width.  If  only  one 
weapon  hit  the  effective . runway  then  the  evaluation  stops 
here.  If  there  are  more  than  one  weapon  hits  on  the 
effective  runway,  then  the  distances  between  two  adjacent 
hit  locations  are  compared  against  the  effective  takeoff 
width.  If  any  of  the  distances  is  greater  than  the 
Ik  effective  takeoff  width,  then  the  runway  is  not  cut.  If 

all  the  distances  are  less  than  the  effective  takeoff  width 

then  the  runway  is  cut  for  this  one  "attack”.  The  "attack" 

\ 

is  repeated  many  times  to  obtain  an  estimate  of  the 
probability  of  cut.  This  results  in  a  binomial  random 
variable  (i.e.,  runway  is  cut  or  not  cut)  and  for  a  large 
number  of  iterations  (attacks"  the  proportion  of  times  the 
runway  is  cut  can  be  approximated  by  a  normal  distribution 
to  get  an  estimate  on  the  number  of  iterations  necessary  to 
achieve  a  desired  accuracy  in  the  estimation  of  the 
probability  of  cut.  To  obtain  a  reasonably  accurate  PC 
estimate  using  the  simulation,  the  number  of  iterations 
(NITERA)  necessary  is  estimated  as  follows: 

Let  p  equal  the  true  probability  of  cut.  If  NITERA 
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equals  the  number  of  trials  and  0  equals  the  number  of 
successful  cuts  observed,  Bernoulli's  theorem  says  that  the 
difference,  d,. between  the  proportion  of  successes  in  NITERA 
trials  and  the  true  probability  of  success  in  a  single  trial 
tends  to  zero,  as  NITERA  approaches  infinity.  Another  way  of 
saying  this  is  the  relation: 


NITERA  P  ^  d 


(15) 


as  NITERA  d-*0. 

We  wish  to  determine,  an  estimate  of  the  true  probability 
of  success  such  that 


where  6/NITERA  is  our  estimate  of  p  and  1 -c£  is  the 
probability  that  our  estimate  does  not  deviate  from  p  by 
more  than  d.  If  NITERA  is  large  enough,  then. the  binomial 
distribution  can  be  approximated  by  a  normal  distribution. 
Using  this  approximation  we  can  show  that: 


(17) 


where  z<^/2  ^ke  two-tailed  standardized  normal  statistic 
for  the  probability  we  seek.  (Ref  13:  191-2) 

In  order  to.  calculate  NITERA  we  need  an  estimate  of 
p  and  q.  If  we  use  the  required  prdbability  of  cut  as  the 
estimate  of  p  and  one  minus  that  as  the  estimate  of  q  then: 
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For  a  99%  confidence  that  the  true  probability  of  cut, 
p,  lies  within  _+  .01  of  the  calculated  probability  PC.  we 
have 


p  =  .95635 
q  =  .04365 
d  =  .01 

z 2  “  2. 58  • 

and  so 

WERA  =  .(2, 58)2( .  95635  ).(.  04365  )  =  2,778  (18) 

( .  01  )2 

With  this  large  number  of  iterations  necessary  each  time  the 
simulation  subroutine  calculates  a  probability,  the 
searching  process  consequently  becomes  very  time  consuming. 
Choice  of  Evaluation  Routine 

There  are  two  evaluation  routines  In  the  model.  The 
first  is  a  discrete  approximation  routine  which  is  very 
efficient  for  the  range  of  input  variables  identified 
earlier.  The  second  routine  is  the  monte-carlo  simulation. 
The  discrete  approximation  is  used  for  problems  where  the 
difference  between  the  effective  runway  width  (ERUNWl)  and 
twice  the  effective  takeoff  width  (ETOWID)  is  less  than 
twice  the  STEP  size.  For  problems  that  do  not  meet  this 
criterion,  the  monte-carlo  simulation  is  used.  This  routine 
is  not  preferred  because  of  the  amount  of  computer  time 
required.  Because  of  the  high  probabilities  of  cut  that  are 


estimated,  around  0.956,  the  number  of  iterations  necessary 
to  be  99^  confident  that  the  actual  probability  of  cut 
lies  within  +.  .01  of  the  estimated  value  requires 
approximately  2,800  iterations.  The  concern  for  speed  will 
by  an  important  factor  in  the  next  chapter,  where  the  dearch 
routine  evaluates  many  different  strategies  in  its  search 
for  the  ’'best”  strategy. 


IV.  Search  for  the  "Best"  Strateg 


The  first  step  in  the  search  is  to  define  the  minimum 
probability  of  cut  for  each  cut  that  will  assure  an  overall 
probability  of  closure  that  meets  or  exceeds  the  required 
probability  of  closure.  This  is  done  by  taking  the  required 
probability  of  closure  and  taking  the  (NCUTS)  root,  where 
NCUTS  represents  the  number  of  cuts  necessary  to  close  the 
runway.  NCUTS  is  the  number  of  shortest  runway  lengths 
(SRL)  in  the  runway  with  a  six  sigma  overlap.  (Refer  to 
chapter  III,  figure  4)  The  derivation  of  the  SRL  is  given 
in  appendix  C.  The  overlap  provides  a  three  sigma  distance 
from  the  end  of  the  SRL  for  each  cut,  in  order  to  have 
independence.  This  minimum  value  for  each  cut  is  designated 
PCSTAR  and  is  the  minimum  requirement  for  the  probability 
of  cut  (PC)  for  each  cut  used  in  the  search  algorithm.  The 
algorithm  searches  for  the  "best"  strategy,  in  terms  of  aim 
points,  locations  of  aim  points  and  number  of  weapons  per 
aim  point,  which  yields  the  highest  PC  for  a  single  cut. 
Results  for  the  entire  runway  are  calculated  using  all  the 
cuts  required  where  the  results  of  each  cut  are  independent 
of  each  other.  For  example,  the  probability  of  closure  for 
NCUTS  identical  cuts  is 

PCL  =  (PC)NGUTS  (19) 

and  the  number  of  weapons  necessary  for  this  strategy  is 

N  =  (n) (NCUTS )  (20) 
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where  is  the  probability  of  runway  closure  using 

identical  cuts.  NCUTS  represents  the  number  of  cuts 
required.  -N  .represents  the  total  number  of  weapons  to 
achieve  the  P^  calculated  above  and  n  represents  the 
number  of  weapons  for  the  ’’best"  strategy  for  each  cut. 

The  problem  is  to  find  the  "best"  strategy  i.e.,  the 
one  which  uses  the  least  number  of  weapons  to  achieve  a 
probability  of  cut  at  least  as  great  as  PCSTAR. 

A  strategy  is  defined  by  three  variables:  NUMAIM, 
NUMBOM(I),  AIMLOC(I)  where 

NUMAIM  -  specifies  the  number  of  aim  points  under 
consideration. 

NUMBOM(I)  -  specifies  the  number  of  weapons  per  aim 
point  for  all  NUMAIM  aim  points. 

AIMLOC(I)  -  specifies  the  location  of  each  aim  point 
for  all  NUMAIM  aim  points,  measured  in  feet  from  one  edge 
of  the  runway. 

The  first  step  is  to  determine  a  lower  bound  and  an 
upper  bound  on  the  number  of  weapons  required.  The 
minimum  number  (MIN)  is  found  by  calculating  the  number  of 
weapons  needed  if  there  were  no  variation  in  weapon  impact 
points  (CEP=0).  This  is  done  by  counting  the  number  of 
non-overlapping  minimum  takeoff  width  sections  within  the 
runway  width.  MIN  is  equal  to  the  greatest  integer  not 
exceeding  ERUNWI/ETOWID.  MIN  represents  a  minimum  number 
of  weapons  necessary  to  cut  the  runway.  When  the  actual 
CEP  is  considered, .  if  MIN  is  feasible  (PC  greater  than  or 


*  :* 


equal  to  PCSTAR),  then  the  least  number  of  weapons  is  MIN; 
if  not,  the  minimum  feasible  soultion  will  not  be  less  than 

MIN. 

A  maximum  number  of  weapons  necessary  to  cut  the  runway 
(MAX)  is  a  feasible  solution  (achieves  a  PC  greater  than  or 
equal  to  PCSTAR. )  From  the  calculation  of  the  MIN  number 
of  weapons  we  also  have  the  minimum  number  of  aim  points 
necessary  to  cut  the  runway.  These  MIN  number  of  aim 
points  are  evenly  spaced  across  the  runway  at  locations 
ERUNWl/ (MIN+1 ) .  This  number  of  aim  points  and  the  aim 
point  locations  are  fixed  in  determining  MAX.  The  number 
of  weapons  per  aim  point  will  be  the  only  factor  varied  to 
achieve  a  feasible  solution.  This  will  be  done  by  first 
allocating  one  weapon  per  aim  point,  evaluating  the  strategy 
allocating  two  weapons  per  aim  point,  evaluating  the 
strategy,  and  so  forth  until  the  PC  is  greater  than  or 
equal  to  PCSTAR.  The  desired  minimum  feasible  number  of 
weapons  will  lie  between  MIN  and  MAX. 

The  second  step  is  to  find  the  highest  PC  for  each 
number  of  weapons,  between  MIN  and  MAX,  and  to  stop  when 
the  PC  calculated  equals  or  exceeds  PCSTAR,  as  shown  in 
figure  12,  where  each  point  represents  the  highest  PC  for 
a  given  number  of  weapons.  With  the  number  of  weapons 
constant,  there  are  many  combinations  of  NUMAIM,  NUMBOM(l), 
and  AIMLOC(I)  that  are  possible.  Because  of  the  symmetric 
nature  of  the  problem,  several  assumptions  will  be  made  j.n 
order  to  reduce  the  number  of  combinations  to  be  considered. 
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Figure  12.  PC  versus  Number  of  Weapons 
Assumption  1  -  Symmetric  numbers  of  weapons  about  the 
runway  center  are  better  than  asymmetric  numbers  of  weapons 
Assumption  2  -  Symmetric  locations  of  aim  points  about 
the  runway  center  are  better  than  asymmetric  locations. 
Using  these  assumptions,  it  is  possible  to  assign  weapons 
to  aim  points  .in  a  systematic  manner,  > Given  the  number  of 
weapons  and  the  number  of  aim  points,  these  combinations 
can  be  easily  identified,  and  are  listed  for  one  to  nine 
weapons  in  table  I.  For  example,  if  the  number  of  aim 
points  is  one,  then  all  the  weapons  will  be  aimed  at  the 
center.  If  the  number  of  aim  points  is  two  then  the  number 
of  weapons  will  be  evenly  divided  (if  possible)  between  the 
two  aim  points.  For  some  cases  there  is  more  than  one 
combination  possible  for  a  given  number  of  weapons  and 
number  of  aim  points.  For  example,  six  weapons  to  be 
distributed  among  three  aim  points  can  be  combined  either 
(1,4*1 )  or  (2,2,2).  Each  of  these  combinations  is  readily 


Table  I.  Combinations  of  Weapons  and  Aim  Points 


Number  of  Aim 

Points 

1 

2  3 

4 

5 

1  J21 

2 

3  (3) 

4  Til 

5  _(5)_ 

6  (6) 

7  (7) 

8  (8) 

9  (9) 


(1,1,1) 


(2,2)  (1,3,3) 


(1,3,1) 


(3,3)  (1,4,1) 

(2,2,2) 


(1,1, 1,1) 


(1,1, 1,1,1) 


(4,4)  (1,6,1) 
(2,4,2) 
(3,2,3) 


(1.7.1) 

(2.5.2) 

(3.3.3) 

(4.1.4) 


(1,2, 2,1)  (1,1, 2, 1,1) 

(2, 1,1, 2) 


(1,1, 3, 1,1  ) 
(1 ,2,1 ,2,1  ) 
(2, 1,1, 1,2) 


(1.3. 3.1)  (1,1, 4, 1,1) 

(2. 2. 2. 2)  (1,2, 2, 2,1) 

(3, 1,1, 3)  (2, 1,2, 1,2) 


(1.1. 5. 1.1) 

(1.2. 3. 2.1) 

(1.3. 1.3.1) 

(2. 1.3.1 .2) 
(2, 2, 1,2, 2) 
(3, 1,1, 1,3) 


identified  and  is  a  possible  combination  for  consideration 
With  the  NUMAIM  and  the  NUMBOM(l)  chosen,  what  are  the 
"best"  locations  for  these  aim  points?  The  aim  points  are 
located  symmetrically  with  respect  to  the  runway  center, 
with  the  odd  aim  point  (if  there  is  one)  being  located  in 
the  center  of  the  runway. 

Because  of  symmetry,  we  only  need  to  specify  half  of 


at 


Center 

1  by  symmetry 

; 

1  aim  point 

Figure  13.  Location  of  Symmetric  Aim  Points 
the  aim  point  locations,  with  the  other  half  being  equal  to 
the  runway  width  minus  the  aim  point  location  of  the 
corresponding  symmetric  aim  point.  For  example,  if  there 
are  two  aim  points,  aim  point  one  is  located  50  feet  from 
the  edge  of  the  runway,  then  aim  point  two  is  located 
(RUNWID-50)  feet  from  the . edge  of  the  runway.  Thus,  the 
two  aim  points  can  be  thought  of  as  being  one  aim  point 
pair.  See  figure  13.  The  search  is  reduced  to  half  of 
the  aim  point  locations  on  half  of  the  runway. 

The  response  curve,  probability  of  cut  versus  distance 
from  the  runway  edges  for  a  pair  of  symmetric  aim  points, 
was  shown  to  be  unimodal  by  Pemberton  (Ref  11:31).  This 
means  that  PC  is  a  non-decreasing  function  up  to  the  "best" 
location,  non-increasing  past  the  "best"  location  and  there 
is  only  one  "best"  location  for  each  pair  of  aim  points 
when  all  other  aim  points  are  fixed.  See  figure  14. 

The  search  for  two  "best"  locations  starts  with  aim 
point  one  at  zero  (one  edge  of  the  runway),  aim  point  two, 
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Figure  14*  Response  Curve  for  One  Pair  of  Aim  Points 
by  symmetry,  at  the  other  edge  of  the  runway  and  evaluating 
the  PC  for  this  strategy.  Next,  the  aim  points  are  moved 
five  feet  toward  the  center  of  the  runway,  evaluated,  and 
if  this  PC  is  higher  than  the  previous  PC  the  movement 
continues  until,  either  the  PC  decreases,  indicating  we 

have  passed  the  "best"  location,  or  the  aim  points  converge 

\ 

in  the  center  of  the  runway,  indicating  one  "best"  location. 
When  the  PC  decreases,  the  aim  points  are  repositioned  back 
to  the  second  to  the  last  aim  point  locations,  and  the 
search  is  repeated  using  one  foot  movements  toward  the 
center  of  the  runway  until  the  PC  again  decreases.  When 
this  occurs  again,  the  aim  point  locations  which  gave  the 
highest  PC  value  is  saved  as  the  "best"  locations  for  two 
aim  points.  The  search  for  three  "best"  aim  point  locations 
is  identical  to  the  search  for  two  "best"  locations,  with 
the  middle  aim  point  fixed  in  the  center  of  the  runway. 

The  search  for  four  "best"  aim  point  locations  starts 


Figure  15.  Initial  Aim  Point  Locations  for  4  Aim  Points 
with  the  four  aim  points  evenly  spaced  across  the  runway 
width.  See  figure  15.  Aim  points  1  and  4  will  be  referred 
to  as  the  outer  pair,  and  aim  points  2  and  3  will  be  the 
inner  pair.  The  search  uses  the  unimodal  idea  for  a  pair 
of  aim  points  by  fixing  one  pair  and  finding  the  "best” 
locations  for  the  other  pair.  Initially,  the  outer  pair 
are  fixed  and  the  "best"  locations  for  the  inner  pair  are 
found.  Next,  the  inner  pair  are  fixed  at  the  "best" 
locations  and  the  "best"  locations  for  the  outer  pair  of 
aim  points  are  found.  Next,  the  outer  pair  is  fixed  at 
its  "best"  locations  and  new  "best"  locations  are  found 
for  the  inner  pair.  This  back  and  forth  procedure  stops 
when  the  "best"  locations  for  both  pairs  of  aim  points  do 
not  change  more  than  one  foot.  This  allows  the  aim  points 
to  "home  in"  on  the  "best"  locations.  The  search  for  five 
"best"  locations  is  identical  to  the  search  for  four  "best" 


locations,  with  the  middle  aim  point  fixed  in  the  center  of 
the  runway.  A  flow  chart  of  the  search  routine  is  presented 


in  figure  16. 

The  search  begins  with  MIN  weapons,  finds  the  "best" 
number  of  aim  points,  "best"  locations  for  these  aim  points 
and  the  highest  PC  for  the  "best"  strategy.  The  "best" 
number  of  aim  points  is  the  smallest,  number  of  aim  points 
that  contribute  to  increase  PC.  If  the  PC  with  one  aim 
point  is  higher  than  the  highest  PC  with  two  aim  points 
then  the  search  stops  for  that  number  of  weapons.  If  the 
highest  PC  with  two  aim  points  is  higher  then  it  is 
compared  to  the  highest  PC  with  three  aim  points,  etc.  The 
"best"  locations  for  pairs  of  aim  points  are  found  as  before. 
Once  the  highest  PC  has  been  found  for  a  given  number  of 
weapons,  the  PC  is  compared  to  PCSTAR.  If  PC  is  less  than 
PCSTAR,  the  number  of  weapons  is  increased  by  one  and  the 
search  begins  again.  If  PC  is  greater  than  or  equal  to 

PCSTAR,  then  the  "best"  strategy  has  been  found  and  the 

\ 

search  stops. 

The  next  chapter  gives  an  overall  description  of  the 
computer  program. 
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Program  Description 


A  computer  program  called  RAM  (Runway  Attack  Model)  was 
developed  to  aid  in  planning  a  runway  attack.  The  diagram 
in  figure  17  shows  the  overall  structure  of  RAM.  Three 
separate  modes  of  operation  are  available  in  the  program. 
Inputs  common  to  all  three  modes  are: 

Runway  dimensions 

-  runway  length  (RUNLEN)  in  feet 

-  runway  width  (RUNWID)  in  feet 
Minimum  launch  window  (MLW) 

-  MLW  length  (TOLEN)  in  feet 

-  MLW  width  (TOWID)  in  feet 
Weapon  characteristics 

-  yield  (W)  in  pounds  TNT 

-  accuracy  (CEP)  in  feet 

where  the  input  units  of  W  in  pounds  TNT  is  converted  to 
the  weapon  damage  radius  by  the  relation: 


R  =  3.54(W) 


(Ref  11 :  58) 


(21) 


and  CEP  in  feet  is  converted  to  sigma  as  derived  in 
appendix  B. 

Mode  1.  In  this  mode  the  user  specifies  the  desired 
probability  of  runway  closure  (PqL»  0<PCIiT).  The  program 
then  searches  for  the  "best”  strategy,  with  identical 
strategies  per  cut,  that  will  give  an  overall  probability 
of  runway  closure  equal  to  or  greater  than  the  desired 
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probability. 

Mode  2.  In  this  mode  the  user  specifies  the  number  of 
weapons  available  (NUMBER).  The  program  then  searches  for 
the  strategy  that  gives  the  highest  probability  of  runway 
closure. 

Mode  3.  In  this  mode  the  user  specifies  an  independent 
strategy  in  terms  of: 

Number  of  aim  points  for  each  cut 
Number  of  weapons  for  each  aim  point 
Location  of  each  aim  point 

with  the  number  and  location  of  cuts  the  same  as  that 
identified  in  mode  1 .  The  program  then  uses  the  appropriate 
evaluation  routine  and  gives  the  approximate  probability  of 
closure  for  the  independent  strategy.  This  mode  does  not 
use  the  search  routine. 

As  an  example  of  this  model,  the  following  were  used 
as  input: 

Runway 

-  length  8,000  feet 

-  width  150  feet 

Minimum  launch  window 

-  length  2,000  feet 

-  width  50  feet 

Weapon  characteristics 

-  yield  250  lbs  TNT 

-  CEP  20  feet 

This  problem  required  a  search  for  a  strategy  to  deny  a 


space  equal  to  2,000  by  50  feet  in  a  runway  of  8,000  by 
150  feet.  Each  weapon  had  a  damage  radius  of  22.3  feet 
and  an  accuracy  of  20  feet  CEP.  For  mode  1,  the  desired 
was  input  as  0.8.  The  resulting  "best"  strategy  was 
to  cut  the  runway  at  four  locations,  1750,  3250,  4750  and 
6250  feet  from  one  end  of  the  runway.  For  each  cut,  three 
aim  points  were  identified  at  36,  75  and  114  feet  from  one 
edge  of  the  runway.  The  number  of  weapons  for  each  aim 
point  was  one.  With  this  "best"  strategy  the  level  of 
runway  closure  was  0.91*  and  the  total  number  of  weapons 
necessary  was  12.  Refer  to  figure  18  for  sample  output. 

In  mode  2,  the  number  of  weapons  available  was  input 
as  15.  Because  of  the  same  input  conditions,  the  number 
of  cuts  will  not  change  for  this  example.  The  number  of 
available  weapons  were  equally  divided  among  the  cuts  with 

the  remaining  weapons  allocated  one  to  a  cut  starting  at 

\ 

one  end  of  the  runway.  For  this  example,  three  cuts  had 
four  weapons  and  one  cut  had  three  weapons.  The  program 
then  determines  the  "best"  strategy  to  employ  for  the 
number  of  weapons  allocated  to  each  cut.  In  this  case  the 
"best"  strategy,  for  cuts  with  four  weapons,  was  to  have 
two  weapons  each  on  two  aim  points  located  at  44  and  106 
feet  from  one  edge  of  the  runway.  For  the  one  cut  with 
three  weapons  the  same  strategy  as  found  in  mode  1  is  the 
"best"  where  one  weapon  was  on  each  of  three  aim  points 
located  at  36,  75  and  114  feet  from  one  edge  of  the  runway. 
The  highest  probability  of  closure  with  15  weapons  availabl 
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Mode  1  output 
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was  0.96.  Refer  to  Figure  19  for  sample  output. 

In  mode  3*  the  user  can  obtain  an  evaluation  of  any 
independently  arrived  at  targeting  strategy.  The  user  may 
want  to  compare  his  strategy  with  the  "best”  strategy  found 

in  either  mode  1  or  mode  2.  The  user  must  use  the  same 

number  of  cuts  and  the  same  cut  locations  but  is  free  to 
specify  the  total  number  of  weapons,  the  number  of  aim 
points  per  cut,  the  number  of  weapons  per  aim  point  and  the 
aim  point  locations.  The  independent  strategy  for  this 
example  was: 

4  Cuts  (as  identified  in  mode  1) 

Cut  1  at  1750  feet  with  2  aim  points 
aim  point  1  at  50  feet  with  1  weapon 
Cut  2  at  3250  feet  with  3  aim  points 
aim  point  1  at  25  feet  with  1  weapon 

aim  point  2  at  75  feet  with  2  weapons 

aim  point  3  at  125  feet  with  1  weapon 
Cut  3  at  4750  feet  with  2  aim  points 

aim  point  1  at  50  feet  with  2  weapons 
aim  point  2  at  100  feet  with  2  weapons 
Cut  4  at  6250  feet  with  2  aim  points 

aim  point  1  at  50  feet  with  1  weapon, 
aim  point  2  at  100  feet  with  1  weapon. 

The  program  evaluates  the  independent  strategy  with  either 
the  discrete  approximation  or  the  monte-carlo  simulation, 
depending  on  the  geometry  of  the  problem  as  discussed  in 
Chapter  III.  The  calculated  probability  of  runway  closure 
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Figure  19.  Mode  2  output 


was  0.68.  Refer  to  figure  20  for  sample  output.  Although 

the  independent  strategy  had  the  same  amount  of  weapons,  12, 

the  differences  in  strategy  produced  a  0.23  difference  in 

the  expected  probability  of  runway  closure. 

Thus,  with  these  three  modes,  a  planner  can  initially 

select  the  desired  level  of  runway  closure.  If  the  number 

of  weapons  required  to  obtain  this  level  is  considered  too 

» 

large,  then  what  level  of  runway  closure  can  be  obtained 
with  a  fixed  amount  of  weapons?  Or,  if  the  "best"  strategy 
is  not  possible  for  a  particular  runway  attack  problem,  due 
to  airfield  defenses  or  terrain,  what  level  of  runway 
closure  can  be  obtained  with  a  modified  strategy?  It  is 
hoped  that  in  answering  these  three  questions,  a  planner 
can  select  the  "best"  strategy  that  meets  the  requirements 
of  a  particular  problem.  For  a  more  detailed  guide  to 

this  program,  refer  to  the  user's  guide  in  appendix  A. 

\ 

The  next  chapter  discusses  the  verification  and 
validation  of  the  model. and  computer  program. 
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VI.  Verification  and  Validation 


Verification  and  validation  are  two  related  processes 
which  increase  confidence  in  the  model.  Verification  is  the 
process  of  making  sure  that  the  model  does  what  the  analyst 
intends  for  it  to  do  and  validation  is  the  process  of  making 
sure  that  the  model  represents  the  real  world.  (Ref  13:208) 

The  verif ication/validation  process  for  Jhe  strategy 
evaluation  phase  used  in  this  research  was  to  develop  a 
reliable  reference  model  (monte-carlo  simulation)  and 
compare  the  new  model  (discrete  approximation)  to  the 
simulation.  The  validation  process  for  the  entire  model 
used  in  this  research  was  to  compare  the  results  from  this 
model  to  the  previous  model  developed  by  Pemberton. 

The  simulation  routine  was  verified  by  running  traces 
for  selected  iterations  and  the  results  were  compared  to 

manual  calculations.  The  simulation  models  the  real  world 

\ 

by  actually  producing  impact  points  along  a  cut  from  the 
underlying  normal  error  distributions.  Then  these  impact 
points  were  ordered  and  the  spaces  between  adjacent  points 
and  the  runway  is  calculated  and  compared  to  the  minimum 
takeoff  width  to  determine  a  runway  cut.  These  impact 
point  spaces  were  calculated  manually  and  compared  to  the 
minimum  takeoff  width.  The  determination  of  a  runway  cut 
agreed  with  the  manual  calculations.  Next,  the  simulation 
routine  was  compared  to  the  discrete  approximation  routine 
and  the  values  from  the  discrete  approximation  were  within 
the  99$  confidence  interval  from  the  simulation  results. 
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See  columns  and  SIM  in  table  II. 

Next,  the  search  routine  was  verified  for  obvious  cases. 
-When  the  minimum  launch  widths  have  large  overlap,  one 
intuitively  expects  that  an  aim  point  at  the  center  of  the 
runway  is  the  "best"  and  this  agrees  with  the  model  where, 
for  a  minimum  launch  window  of  100  feet  within  a  runway 
width  of  150  feet  the  "best"  strategy  had  one  aim  point  in 
the  center  of  the  runway.  Also,  previous  results  by 
Pemberton  were  available  to  compare  against  this  new  program 
See  table  II  for  comparison  between  Pemberton's  model  and 
RAM  (Runway  Attack  Model)  developed  in  this  research. 
Although  agreement  is  not  100$,  it  is  felt  that  improvements 
in  the  calculation  of  the  shortest  runway  length  for  the 
MLW,  the  increased  efficiency  of  the  discrete  probability 
calculation  and  the  new  search  algorithm  are  responsible 
for  the  minor  differences  in  the  number  of  weapons  necessary 
to  achieve  the  predetermined  level  of  Runway  closure  (0.8) 
and  the  estimated  probability  of  runway  closure. 

The  validation  process  discussed  above  is  for  the  mode  1 
operation  in  the  program  where  the  program  searches  for  the 
"best"  strategy  to  obtain  the  predetermined  P^  value. 

Modes  2  and  3,  search  for  highest  PC  for  fixed  number  of 
weapons  and  evaluation  of  any  strategy,  are  subsets  of 
mode  1  and  are  felt  to  be  validated  in  the  previous  analysis 


Table  II.  Validation  Results 
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VII.  Sensitivity  Analysis ' 

The  task  of  sensitivity  analysis  is  very  much  an  art. 
rather  than  a  science.  The  sensitivity  of  this  model  to 
a  change  in  one  or  more  input  variables  depends  on  the 
actual  values  of  these  variables  and  on  how  big  a  change 
is  proposed  on  selected  variables.  Since  these  factors  are 
different  for  each  problem,  no  automatic  sensitivity 
routine  was  included  in  the  model.  A  sensitivity  analysis 
can  be  performed  by  changing  the  input  variables  and  re¬ 
doing  the  problem.  An  example  of  the  sensitivity  of  this 
model  for  one  specific  problem  was  performed  by  changing 
the  input  value  of  one  factor  at  a  time.  The  total  number 
of  weapons  was  chosen  as  the  response  variable.  The  first 
step  was  to  run  the  program  with  the  original  input 
variables . 

Input:  Output: 

RUNLEN  -  8000  N  =  12 

RUNWID  -  150 

TOLEN  -  2000 

T0.WID  -  50 

W  -  250 

CEP  -  20 

PCL  -  0.8 

The  second  step  was  to  vary  one  input  variable  at  a  time 
to  find  the  range  over  which  the  response  variable  (total 
number  of  weapons)  did  not  change,  i.e.,  the  response 
variable  is  insensitive  to  changes  within  this  range  for 


VIII.  Conclusions  and  Recommendations 


This  research  effort  resulted  in  a  computer  program 
that  can  aid  a  commander  or  planner  with  the  decision  of  how 
to  "best”  attack  a  runway.  This  program  is  an  extension  of 
the  previous  work  by  Pemberton  (Ref  11).  Improvements  have 
been  made  in  the  efficiency  of  the  computer  program  and  an 
overall  program  was  developed  to  offer  the  user  three  modes 
of*  operation. 

Conclusions 

Thi s  is  the  only  runway  attack  program  that  does  not 
require  the  user  to  have  an  initial  attack  strategv,  but 
rather  this  program  returns  an  attack  strategy  that  is 
considered  "best"  to  close  the  runway  at  a  given  level. 

This  program  allows  a  planner  to  see  what  one  "best"  strategy 
looks  like.  If  the  "best"  strategy  found  in  this  program 
is  not  feasible  due  to  limited  number  of  weapons  or  tactical 
considerations,  then,  using  inodes  2  and  3  iteratively  will 
lead  to  a  "best"  attack  strategy  for  a  particular  problem 
and  estimate  its  probability  of  runway  closure.  Mode  2 
allows  the  user  to  specify  a  fixed  number  of  weapons  avail¬ 
able  and  returns  the  highest  level  of  runway  closure  using 
the  fixed  number  of  weapons  in  a  "best"  way.  Mode  3  allows 
the  user  to  change  the  individual  cut  strategies,  while 
leaving  the  number  and  locations  of  cuts  the  same  as  that 
identified  from  mode  1 .  to  design  a  tactically  feasible 
attack  strategy.  The  program  will  evaluate  this  strategy 
and  return  the  estimated  level  of  runway  closure  for  this 
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independent  strategy. 

By  using  this  program  in  a  series  of  defining  and  re¬ 
defining  the  runway  attack  problem,  the  user  will  gain 
insights  and  experience  about  the  problem  of  attacking 
runways. 

Recommendations 

While  this  research  effort  answered  many  of  the 
questions  posed  in  the  beginning  of  this  report,  it  raised 
and  left  unanswered  other  questions.  Further  work  and 
research  should  be  directed  in  these  areas: 

1.  Adjust  the  computer  program  to  address  two  other 
problem  scenarios: 

(a)  R&D  scenario,  where  the  weapon  yield  and 
accuracy  are  the  main  factors  in  determining  the  "best" 
strategy. 

(b)  Runway  planning  scenario,  where  the  vulner¬ 
ability  of  a  current  or  proposed  base  to  enemy  attack  can  be 
evaluated  with  the  runway  dimensions  as  the  main  factors 
determining  the  "best"  strategy.  The  objective  here  is  to 
minimize  the  highest  probability  of  runway  closure,  subject 
to  cost  constraints. 

2.  Adapt  this  program  to  an  interactive  mode.  The 
program  is  small  enough  to  warrant  use  in  the  interactive, 
rather  than  batch  mode. 

3.  Pursue  the  order  statistics  approach  for  the 
probability  of  cut  presented  in  Chapter  III. 

4.  Expand  the  discrete  approximation  routine  to 


include  cases  where  s^  and  s|^ sp W  2  n0^  overlap. 

5.  Consider  non-identical  cuts  in  obtaining  the 
"best”  strategy  to  obtain  a  certain  level  of  runway  closure 

6.  Add  an  automatic  sensitivity  analysis  routine. 

(a)  Sensitivity  can  be  in  terms  of  changes  to 
total  number  of  weapons  as  in  Chapter  VII. 

(b)  Sensitivity  analysis  taking  two  factor  and 
higher  interactions,  into  consideration. 

(c)  Analyzing  the  sensitivity  of  the  "best" 
strategy  to  changes,  in  input  variables. 

7.  Consider . non-circular  error  distributions  for 
weapons  in  terms  of  REP  and  DEP. 
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User's  Guide  to  Runway  Attack  Model 


This  is  an  expected  value  model  that  gives  the  "best” 
strategy  to  attack  a  runway.  This  program  is  written  in 
Fortran  77,  requires  approximately  63,300  bytes  of  central 
memory  ( CDC  6600)  and  takes  approximately  10  seconds  of 
execution  time  to  run  an  average  problem.  The  execution 
time  requirement  increases  to  approximately  200  seconds 
when  the  simulation  approximation  subroutine  is  being  used. 
The  program  uses  calls  to  the  International  Mathematical 
and  Statistics  Libraries  (IMSL)  for  procedures: 

MDN0R  -  returns  the  cumulative  normal  distribution 
GGNML  -  returns  a  normal  random  variate 
VSRTA  -  returns  a  sorted  array. 

This  program  has  three  modes. 

Mode  1 

In  this  mode  the  user  specifies  the  desired  level  of 
runway  closure  and  the  program  returns  the  "best"  strategy 
to  achieve  at  least  that  level  of  runway  closure. 


Mode  2 

In  this  mode  the  user  specifies  the  number  of  weapons 
available  and  the  program  returns  the  "best"  strategy 
to  use  and  gives  the  expected  level  of  runway  closure. 


In  this  mode  the  user  specifies  an  independent  strategy, 
using  the  number  and  location  of  cuts  identified  from  mode  1 
and  the  program  returns  the  expected  level  of  runway  closure 


Input  variables  are: 


Mode 

1  -.User  specifies  desired  probability 

2  -  User  specifies  number  of  weapons  available 

3  -  User  specifies  attack  strategy 
RUNLEN  -  Runway  length  in  feet 
RUNWID  -  Runway  width  in  feet 

TOLEN  -  Minimum  takeoff  length  in'  feet 
TOWID  -  Minimum  takeoff  width  in  feet 
W  -  Weapon  yield  in  lbs.  TNT 
CEP  -  Weapon  accuracy  in  feet 

Pd  “  Mode  1 ,  desired  probability  of  runway  closure 
NUMBER  -  Mode  2,  number  of  weapons  available. 

Card  Deck  Setup  for  the  CDC  6600 


job  card 


FTN5. 

ATTACH, IMSL, ID=LIBRARY, SN=ASD. 
LIBRARY, IMSL. 

LGO. 

*E0R 


Program 


*E0R 

*E0R 


Data 


6  Q 


Data  Cards: 

Card  1  mode  (1 ,  2,  3) 

Card  2  RUNLEN,  RUNWID,  TOLEN,  TOWID,  W,  CEP 

Card  3  PCL  or  NUMBER 

Card  3  and  on  (mode  3  operation  specifying  an  independent 
strategy) 

-  number  of  cuts  (as  identified  from  mode  1) 

-  location  of  cut  1 ,  number  of  aim  points  for 
cut  1 

-  aim  point  location  1,  number  of  weapons  for 
aim  point  1 

-  aim  point  location  2,  number  of  weapons  for 
aim  point  2 


-  location  of  cut  2,  number  of  aim  points  for 
cut  2 

-  aim  point  location  1 ,  number  of  weapons  for 

aim  point  1 

-  aim  point  location  2,  number  of  weapons  for 

aim  point  2 


-  location  of  last  cut,  number  of  aim  points 

for  last  cut 

-  aim  point  location  1,  number  of  weapons  for 


aim  point  location  1 


-  aim  point  location  2,  number  of  weapons  for 
aim  point  location  2. 

Mode  1  example: 

1 

8000.  150.  2000.  50.  1000.  50. 

0.8 

*E0R 

This  card  set  up  will  return  the  "best"  strategy  for  closing 
a  runway  (8,000  by  150  feet)  for  a  MLW  (2,000  by  50  feet) 
with  weapons  (1,000  lbs  TNT,  50  foot  CEP)  with  a  desired 
level  of  closure  0.8.  See  figure  A-1  for  results. 

Mode  2  example: 

2 

8000.  150.  2000.  50.  1000.  50. 

15 

*E0R 

This  card  set  up  will  return  the  "best"  strategy  for 
employing  15  weapons  against  the  same  runway  as  above. 

See  figure  A-2  for  results. 

Mode  3  example: 

3 

8000.  150.  2000.  50.  1000.  50. 


This  card  set  up  will  return  the  expected  probability  of 
runway  closure  for  this  independent  strategy.  Note,  the 
number  of  cuts  used  and  cut  locations  are  identical  to  the 
output  from  mode  1  as  shown  in  figure  A-1 .  Note  also, 
there  is  no  requirement  for  symmetry  for  the  aim  point 
locations,  or  number  of  weapons  per  aim  point  for  this  mode 
of  operation.  See  figure  A-3  for  results. 
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Figure  A-1 .  Mode  1  output 


RUEWAY 

fIMMliM  LAUECH 

WEAPON 

NUVRER 

u:ed:w 

r  F 

LENGTH  WIDTH 

LENGTH  WIDTH 

YIELD 

CEP  WEAPON 

5CCC.  150. 

2  G  0  C •  5  0. 

1 G  C  0  . 

50.  15 

ELir-  WAY 

0 f. 0  •  BY 

150 

•  FEET 

VI f  LAUNCH  WT 

EDOW  rote.  BY 

50 

.  FEET 

WEAPON  CHARACTERISTIC 

YIELD  lCLO.iiC  FGUNC' 
CEP  50 •  FEET 

P'CBABILITY  CF  CL33LFE  .76 

TOTAL  NUKBEF  CF  WEAPON  15 


AIN  POINT 

3 

NL'KBc  ; 

(LENGTH* 

WIDTH* 

WEAPCNS 

1 7  5  0  •  iw  . 

54.0  5 

2 

17«0  .0C 

6.0  C 

c. 

3250. Ct 

C*.CC 

2 

3250.CL 

'  6  .C  C 

£ 

4  750  .c : 

-4.JC 

4  750 •  )  I 

6. CO 

2 

6  2 1: 0  »  0  C 

56.0C 

1 

6250 .CC 

75.CC 

1 

62 5C.CT 

4.C3 

1 

Figure  A-2.  Mode  2  output 


RUNWAY  MIfilVUM  LAL'jCH  WEAPON 

Wlf.D.n 


LET GTH 

WIDTH 

LENGTH  W^OTH 

YIELD 

CEP 

PLOD. 

15  0. 

20 OC.  EC. 

1 C  0  C  . 

5C  • 

INDEPENDENT  ANALYSIS 
INDEPENDENT  strategy: 

CUT  AIK  POINT  .NtNBEP  PE' 

LC  CAT  I CN  LOCA'ICN  AIM  PCIf.'r 


75.00  4 

3  2  5  C  • 

50. CO  3 

75. OC  ? 


25. CO  1 

75. CO  2 

100.  C  O  2 

625L. 

25.00  1 

5C.CC 

75.00  2 

i  o  c . :  o  l 

PROBABILITY  OF  RUNWAY  CLOSURE  IS  .90 
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APPENDIX  B 


Relationship  Between  CEP  and  Sigma 
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Relationship  Between  CEP  and  Sigma 

The  weapon  error  distribution  in  this  research  is 
assumed  to  be  circular  normal,  centered  at  the  aim  point. 
The  circular  normal  probability  density  function  is: 


f(x,y)  = 


1 


271a 


exp 


•-i  ~ 

* 


2±  2 
Jx_ 
0-2 


( B-1  ) 


(Ref  11 :  60) 

The  probability  that  the  weapon  will  land  in  a  circle  with 
radius  R  and  area  A  is  given  by: 


P(A)  = 


1 


il 


f^)j 


dx  dy 


2!T<r  A  exPl 
Transfering  to  polar  coordinates  by  letting 

x  =  r  cos  0 
y  =  r  sin  0 


(B-2) 


yields 


P(A)  "ij  ■{{-{•‘fell' dr  de 

=  J  h[h  l  d9J r  exp{-r?/s°2  ]dr 

=  j  r/<y2  exp{-r2/2<y2jdr 

)J-R2/2<r2} 


(B-3) 


=  1  -  exp< 

i 

CEP  is  defined  as  the  radius  of  a  circle  with  a  0.5 
probability  of  a  weapon  landing  within  it,  or: 


P ( A )  =  0.5  =  1  -  exp 


-{cep2/ 2  a2] 


(B-4) 
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moving  terms  from  one  side  to  the  other  yields 


1  -  0.5  =  exp|-CEP2/  2<r] 
0.5  =  expj-CEP2/2o-2} 
ln(0.5)  =  -CEP2/ 2  o'2 


(B-5 ) 


This  relationship  allows  converting  the  normal  units 
of  weapon  accuracy,  CEP,  into  sigma,  which  is  used  in  the 
research. 


Development  of  Shortest  Runway  Length 


The  determination  of  the  minimum  number  of  cuts 
required  to  close  a  runway  involves  finding  the  shortest 
runway  length  that  completely  contains  a  minimum  launch 
window  (MLW,  takeoff  length  by  takeoff  width). 

When  the  weapon  damage  radius  is  taken  into 
consideration,  the  dimensions  for  the  runway  and  the  MLW  are 
increased  by  an  amount  equal  to  the  damage  radius  all 
around.  This  presents  a  problem  because  not  the  shapes  are 
not  rectangular,  rather,  they  have  corners  that  are 
described  by  an  arc  with  radius  equal  to  the  damage  radius. 

In  order  to  find  the  shortest  runway  length  that 
sompletely  contains  a  MLW,  the  MLW  was  approximated  by 
a  rectangle  of  dimensions  effective  takeoff  width  by 
effective  takeoff  length.  The  shortest  distance  was 
approximated  by  the  base  of  the  triangle  formed  by  the 
effective  takeoff  length  and  the  line  prependicular  to 
the  runway  that  touches  the  corner  of  the  MLW.  This 
distance  was  chosen  because,  it  under-estimates  the  actual 
distance.  The  distance  from  one  prependicular  to  the 
next  prependicular  overestimates  the  actual  distance. 

The  underestimation  is  preffered  because  it  is  better  to 
have  too  many  cuts  rather  than  too  few. 

This  estimate,  called  the  shortest  runway  length 
(SRL)  can  be  found  by  using  a  geometric  argument.  From 
figure  B-1,  RW  can  be  divided  into  two  segments  A  and  B. 


RW  -  Effective  Runway  Width 
TL  -  Effective  Takeoff  Length 
TW  -  Effective  Takeoff  Width 
SRL  -  Shortest  Runway  Length 


SRL 


Figure  C-1 .  Shortest  Runway  Length 


The  triangle  A,  SRL,  TL  is  similar  to  the  triangle  C,  TW,  B. 
Since  these  two  triangles  are  similar,  the  ratios 


TL/SRL  =  TW/B 


( B—  1  ) 


are  equal.  Another  relationship  is 


(TL)2  =  (SRL)2  +  (RW-B)2 


(B-2) 


We  now  have  two  equations  in  two  unknowns,  SRL  and  B. 
Solving  for  B  from  equation  ( B— 1 )  yields 


B  =  SRL(TW/TL) 


(B-3) 


Substituting  equation  (B-3)  into  equation  (B-2)  yields 

(TL)2  =  (SRL)2  +  (RW-SRL(TW/TL))2  (B-4) 


snuaring  the  terms  yields 


(TL)2  =  (SRL)2  +  (RW)2  -  2(SRL)N(RW)(TW)/(TL)  + 

(srl)2(tw)2/(tl)2 


(B-5) 


collecting  terms  yields 


[l  +  (TW)2/(TL)2J(SRL)2 
[(RW)2  -  (TL)2]  =  0 


-  [2(RW)(TW)/(TL)J  SRL 


(B-6) 


applying  the  quadratic  formula  yields 


[2 ( RW )  ( TW ) / Tl]  +  V^2(RW)(TW)/Tl]2-U)f  1 +TW2  rv;2_tl2 

L  TL2JL 


2(1 +TW4/TL  ) 


2  2 
TL  +  RW 

=  (rw)2(tl)  z.  tl2Vti? 

2  2 

TL  +  RW 

=  (RW)2(TL)  -  TL3 

2  2 

TL  +  RW 

taking  the  positive  root, 

=  (RW)2(TL)  +  TL3 

2  2 

TL*  +  RW* 

=  TlFtl2  +  RW2] 

2  2 

TL*  +  RW* 

=  TL 

The  solution  is  real  as  long  as  the  determinant 

2  2  2 

TL  +  TW*  -  RW* 

is  positive.  Geometrically,  this  is  saying  that  this 
procedure  works  for  cases  where  the  diagonal  of  the 
minimum  launch  window  is  greater  than  the  effective  runway 
width.  If  this  were  not  so,  the  MLW  could  be  positioned 
vertically  on  the  runway. 
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APPENDIX  D 


PROGRAM  RAM 

DIMENSION  NUMBOM<20 ). AIMLOC(20) 

DIMENSION  NPERCT (  20  ) 

COMMON  PROSEC (20*20 *2 ) 

COMMON  /INPUT/  RUNLEN * RUNUI 0 * TOL E N * TOU I D . U»CEP 

COMMON  /TRANS/  ERUNLE , E RUN WI .ETOLEN . ETOUI 0 .DAHR AD. ST ADEV 

COMMON  /GRP/  NSPU.STEP 

COMMON  /RUAY/  SRL  »PCSTAR»NCUTS»CUTL0C(20> 

COMMON  /BEST/  PCBEST ( 20  ) »NALBST (2 0 ) .MBNBST ( 20 .5 ) « AIMBST (20 *5 > 

COMMON  /SIMS/  DSEED.NITERA 

COMMON  SIMCAL 

DOUBLE  PRECISION  OSEED 

LOGICAL  SIMCAL 

READ* .MODE 

READ* .RUNLEN* RUNUI O.TOLEN.TOWID.M. CEP 

CALL  INIT 

CALL  CONVRT 

IF  (HOOE.EQ.1J  THEN 
READ* .PCLOSE 
PRINT  100 
PRINT  110 
PRINT  120 

PRINT  130 «  RUNLEN «RUNU ID. TOLEN.T QUID. U. CEP. PCLOSE 
CALL  RUNWAY  (PCLOSE) 

CALL  BOUNDS  (MIN  «MAX) 

CALL  SEARCH  (M I N «M A X « NUMB .PC ) 

DO  10  I=1»NCUTS 

NPERCTd  >=NUMB 
CONTINUE 

CALL  EVAL  (  NPEP.CT  .  NUMBER  ,PCLE  ST  > 

CALL  RESLTS  (NPERCT .NUMBER .PCLEST » 

ST  CP 

ELSE  IF  (MOOE.EQ.2)  THEN 
REAO* .NUMBER 
PRINT  102 
PRINT  110 
PRINT  122 

PRINT  1 3 2. R UNLEN. RUNUI D.TOLEN.TOUID.V. CEP* NUMBER 
CALL  RUNUAY  (1  .0  » 

CALL  NUMSET  (NPERCT .NUMBER .PCLEST > 

CALL  RESLTS  (NPEPCT. NUMBER. PCLEST) 

ST  CP 

ELSE  IF  ( HODE.EQ.3)  THEN 
PRINT  103 
PRINT  113 
PRINT  123 

PRINT  130. RUNLEN. RUNUI D» TOLEN.TOUID.W. CEP 
CALL  RUNUAY  (1 .0 ) 

CALL  TNOEP 


STOP 

C 

END  IF 
C 

PRINT  140 
C 
C 
C 

100  FGRMAT!*1*,3X**RUNWAY**5X**MINIMUM  LAUNCH* • 
Z3X**WEAPON**4X«*PROBABILITY*» 

102  F0RMAT!*1*,3X.*RUNWAY*,5X* 'MINIMUM  LAUNCH** 

23X*  *11  CAPON*  #6X*  ‘NUMBER  •  > 

103  FCRHAT!*1*»3X»  *RUNUAY*  *5X  »  *MINIHUM  LAUNCH** 

Z3X*  *UEAPON* ) 

110  FORMAT!  •  •  *  T1  <5  *  •  WINDOW*  *  T47  *  *OF*  I 
113  FORMAT!*  '  * T1 9 « ' WINDOW •  ) 

120  FORMAT! •  LENGTH  UIDTH  LEN6TH  WIDTH*, 4X* 

Z*  YIEL C  CEP • *2  X  *  *CLO  SURE* ) 

122  FORMAT!*  LENGTH  WIDTH  LENGTH  UI0TH*«4X* 

Z* YIELO  CEP*»2X** WEAPONS* ) 

123  FORMAT!*  LENGTH  WIDTH  LENGTH  WIDTH**4X* 

Z* YIELD  CEP • > 

130  FORMAT! *0 *,F6.0,F6.6,2X*F6.0,F6.0 *4X,F5.0 *3X*F4.0 ,3X,F4.2> 
132  FORMA T!*0*«F6.0*F6.0*2X*F6.0*F6.0*4X«F5.0*3X*F4.0*3X ,14 » 
140  FORMAT ! *0 •**  INPUT  ERROR  -  CHECK  MODE  SPECIFIED  »**•> 
END 

SUBROUTINE  INZT 
COMMON  /SIMS/  OSEEO  «NI TERA 
DOUBLE  PRECISION  OSEEO 
C 

OSEEO -58528 SI  24* 

NITER A=2770 

RETURN 

END 

SUBROUTINE  CONVRT 

COMMON  /INPUT/  RUNLEN*RUNWlD«TOLEN*TOWID*W*CEP 
COMMON  /TRANS/  ERUNLE , ERUNWI ,ETOLEN« ETOWI D* OAMR AO • STAOEV 
COMMON  /GRP/  NSPW, STEP 
COMMON  SIMCAL 
LOGICAL  SIMCAL 
C 

0AMRA0=3.54*( (W )» * ! 1./3. > ) 

STAOEV- !CEP)/1 *1774 
C 

STEP=5.0 

ERUNLE=RUNLEN*2.*0APRAD 
ERUNWI=RUNWI0+2.*DAHR AO 
ETOLEN=TOLEN  ♦2.*0AMPAD 
ETOWI D=TOUI D  ♦2.*DAPRA0 
SIMCAL=. FALSE. 

IF!!ERUNWI-2*ET0WID).GT.!2*STEPn  THEN 
SI MCAL-.TRUE. 
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PRINT  100 
PRINT  110 
PRINT  120 
END  IP 

NSPU= C <ERUN VI -ETO VI D>/ STEP >♦. 0001 

FORMAT  <  *0**»  SIMULATION  APPROXHATION  IN  USE  ***•> 

FORMA T<  *0  *• *  RUN  TIME  MAY  NEED  TO  BE  INCREASED  ***•> 

FORMAT <  *0 RESULTS  ARE  ONLY  ACCURATE  TO  ♦  /-  .01  ***•> 

RETURN 

END 

SUBROUTINE  PRONE ( NUMAI M«NUMBOM*AIMLCC .PC) 

DIMENSION  NUMBOM(20 ) « AIMLOCI20 ) «0AML0C(2Q ) 

COMMON  /TRANS/  ERUNLE .ERUNVI »ET0LEN» ETOUI 0 .DAMR AO, STADEV 
COMMON  /ORP/  NSPWt STEP 

DAMLO C( I) =AIMLOC ( 1 )*DAMR AD 
BR= CETONI O-DAMLOC ( 1 > I /ST ADEV 
BL=<ERUNVl-ETOVIO-DAMLOC 111) /STADEV 
IF< BL.OT.BR)  ThEN 

PRINT*#* SI  AND  SN  DO  NOT  OVERLAP. ..DOUBLE  CHECK  INPUTS* 
PC-0 .0 
RETURN 
END  IF 

CALL  MONOR<BR#BONOP> 

CALL  MONOR < 8L« BONOL) 

PC=BONDR-BONOL 

RETURN 

END 

SUBROUTINE  PROCAL  < NUMAIM ,AIMLOC> 

DIMENSION  AIMLOC(20>«DAMLOC(20> 

COMMON  PROSEC <20  #20  #2 ) 

COMMON  /TRANS/  ERUNLE .ERUNWI »ETOLEN»ETOUIO# DAMR AO, STADEV 
COMMON  /GRP/  NSPU.STEP 

DO  10  J=l, NUMAIM 

DAMLOC < J > =AI MLOC < J> +DAMR AD 
CONTI NUE 

00  SO  11=1 #NURAIM 

00  40  12=1 «NSPU 
X=I2-i 

BL=CX*STEP-DAML0C(I1> l/STAOEV 
BR=<X*STEP*ETOVlD-DAMLOC< II >> /STADEV 
CALL  M0N0R<BL«PRL) 

CALL  MONOR(BRtPRR) 

PRR=1.-PRR 
PROS £C<I1#I2*1 )=PRL 
PRO SEC ( I1#I2#2)=PRR 
CONTINUE 
CONTINUE 
RETURN 


SUBROUTINE  UMON(NUHAIH.NUMBOM. AIMLOC.PC) 

DIMENSION  NUM80MI20 ),AIMLCC(20) 

DIMENSION  S(2y20  )  t  Sl(2>,  S2(6> 

COMMON  PROSEC (20 *20  .2 ) 

COMMON  /TRANS/  ERUNLE .ERUNkI .ETOLEN.ETOUIO ,OAMR AD» ST AOEV 
COMMON  /6RP/  NSPU,  STEP 
C 

CALL  PROCAL(NUMAIM.AI MLOC) 

TOTAL  =0  .0 
NH1=NSPU-1 

GAP- (ERUNUI-2*ETOUIO) 

IFCGAP.LE.0.0 )  GO  TO  SO 
IF(GAP«LE.STEP)  GO  TO  30 
IF(GAP.L£.(2*STEP) »  GO  TO  10 
PRINT*.********  USE  SIMULATION  ******** 

RETURN 

10  CONTINUE 

DO  15  1:1.6 

S2(I >  =  1.0 
15  CONTINUE 

DO  20  I1=1.NUMA1M 

SL1=PR0S£C (II. 1.1) 

SL2=PR0$EC(I1 ,2,1 > 

SR NM1= PROSEC (I1.NM1.2) 

SRN= PROSCC < I 1 . NSpy » 2 1 

SM1NH1=PR0SEC( Il.NMl.l)  -  C 1 ,-PR0SEC( I 1 .1 .2 )) 
SM2N=PROSEC< 1 1  .NSPU.l)  -  (1« -PROSEC (I 1.2.2)) 

SM2NM1 =PROSEC( II . NM1 .  1 )  -  ( 1 .-PROSEC ( I 1 ,2 .2 ) ) 

IF  (SM2NM1.LE.  0.)  SM2NM1=0. 

SM3N=PROSEC( 1 1  .NSPU.l)  -  C 1 . -PROSEC ( 1 1 .3.2)) 

IFCSM3N  .LE.  0.)  SM3N=0. 

S2 (1 )  s  (CSL1  ♦  SP1NM1  ♦  SRNM1)**NUKB0H( I1))*S2(1 ) 

S2(2 )  =  ( (SL2  *  SM2N  *  SRN>  ••NUNBOHCI l ) > *S2 (2 > 

S2(3)  =  ( ( SL 1  ♦  SM2NM1  ♦  SRNM1I**NUMB0M( I1))*S2(3) 

S214)  :  ( (SL2  ♦  SM3N  ♦  SRN)  • *NUMBOMC I 1 )) *S2«4 I 

S2(5»  =  ( (SL1  *  SM1NM1  ♦  SRN)  • *NUNBOM( I 1 > ) *S2 (5 ) 

S2 (6 )  s  MSL1  ♦  SH2NM1  ♦  SRN)  •  *NUMBOM(1 1 )  )  *S2(6  ) 

20  CONTINUE 

TOTAL =TCT  AL-S2 (1 )-S2( 2 )*S2 (3)*S2( 4 )*S2 ( 5) -S2(6  > 

30  CONTINUE 

S1<1)=1.0 
Sl(2)=1.0 
DO  40  Ilsl.NUMAIM 

SL=PRGS£C( 11 .1.1 ) 

SM1=PR0SEC(I1 »  NSPU.l)  -  (1  • -PROSEC  ( 1 1  .1.2  )) 
SM2:PR0SEC(I1. NSPU.l)  -  d .-PROSEC U 1 ,2,2 > > 
SR=PR0SEC(I1,NSPU.2) 

IFCSM2  .LE.  0.)  SM2=0. 

Sill)  =  <  (SL  ♦  SHI  ♦  SR)  ••  NUMBOMd  1 ) )  «  Sl(l) 
SK2)  =  ((SL  ♦  SM2  ♦  SR)  •*  NUMBOMd  I ) )  •  SI  (2) 

40  CONTI  NUE 

TOTALsTOT  AL-S1 <1 ) *S 1 (2) 

SO  CONTINUE 
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00  55  1=1 «NSPU 
S(ltl)=1.0 
S(2«I>=1.0 
55  CONTINUE 

DC  70  I2=1«NSPU 

00  60  11=1 tNUM AI M 

S(1»I21  =  <  (PROSECI II  »I  2  » 1  l+PROSECI II *12  »2 >  )  *  •NUHBOM  (1 1 )  )* 
1S(1«I 2) 

IF (I2.EQ.NSPU)  GO  TO  60 
13=12+1 

SC2»I2)=( (PROSECI II ,12,1 )+PROSEC( II »I 3 *2 I I++NUMBOM CIl I > * 
1S(2«I2) 

60  CONTINUE 

70  CONTINUE 

00  80  12=1 «NSPU 

T0TAL=TCTAL+$(1#I2» 

IF(I2.EQ.NSFU)  60  TO  80 
13=12+1 

TOT AL=T0T AL-S ( 2 •! 2  > 

80  CONTINUE 

PC=1. -TOTAL 

RETURN 

END 

SUBROUTINE  SIM<NUMAIM*NUMBOH.AIMLOC*PC) 

DIMENSION  NUMBOMC20 I .AIML0CC2C I tOAMLCC (20 > 

DIMENSION  HIT 1 50) tRl ( 1 ) 

COMMON  /TRANS/  ER UNLC tERUN NI «ETOLEN(ETOUI D .OAMR AOtSTAOEV 
COMMON  /SIMS/  OSEEOtNITERA 
DOUBLE  PRECISION  OSEEO 
C 

Of  EN=O.Q 
INT=1 

DO  10  I =1 vNUMAIM 

OARLOCII  >=AIMLOC(  D+OAMRAO 
10  CONTINUE 

DO  200  ITER=1»NITERA 
NURB=0 

DO  20  J1=1»NUHAIM 
N1=NUMB0M ( J1 ) 

00  20  J2=1«N1 

NUMBsNUMB+1 

CALL  GGNML (DSEED  »INT  «R1 > 

HIT  tNUHBl-OANLOCIUl >  +  (ST  ADEV+R1  (INT  >) 

20  CONTINUE 

CALL  VSRTAIHIT fNUMB) 

NSHOR  T=0 
MLONG=0 
NSHL0=0 

00  30  L=1*NUMB 

IF(HIT(L).LT.O • )  NSHOR T=NSH0RT+1 
IF! HI T ( L ) *6 T*ERUNUI I  NLON6=NLONG+1 
CONTINUE 

IFCNSHORT.EQ.NUHB)  GO  TO  100 
IFCNLONG.EQ.NUMB)  GO  TO  100 
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NSHLO=NSHORT*NLONG 
IFtNSHLO.EQ.MJMB>  GO  TO  100 
C 

C  IS  THE  FIRST  WEAPON  CLOSE  TO  THE  RUNWAY  EOGE? 

C 

NFIRST=NSHORT*l 

IF(HITtNFIRST).GT.ETOWIO)  GO  TO  100 
C 

C  IS  THE  LAST  WEAPON  CLOSE  TO  THE  RUNWAY  EDGE? 

C 

NLAST=NUMB-NLONG 

IFCH!T<NLAST).LT.(£RUNWl-ETOWID>>  GO  TO  100 
C 

C  IS  THIS  THE  ONLY  WEAPON  ON  THE  RUNWAY 

C 

IF(NFIRST.EQ.NLAST)  GO  TO  200 
C 

C  CHECK  ON  THE  OISTANCE  BETWEEN  ADJACENT 

C  IMPACT  POINTS  ON  THE  RUNWAY. 

C 

NLH1=NLAST-1 

DO  40  J=NFIRST .NLMl 

DIS=HIT<J*1 >-HIT(J> 

IFIOIS.6T .ETOWXO)  GO  TO  100 
40  CONTINUE 

C 

C  NO  OPEN  SPACE  FOUN0 

C 

GO  TO  20  0 

100  OPEN=CPEN*l.C 
200  CONTINUE 

XITER  =NIT  ERA 
TOTAL=OPEN/XITER 
PC=  1. -TOTAL 
RETURN 
END  ' 

SUBROUTINE  RUNWAY tPCLOSC > 

COMMON  /INPUT/  RUMEN  »RUNUID«TOLEN*  TOUlOt  W  tCEP 

COMMON  /TRANS/  ERUKLE *ERUNWI .ET0LEN.ET0U1D.DAMR ADtSTADEV 

COMMON  /RWAY/  SRL «PCST AR vNCUTS* CUTL0CC20 > 

C 

SRL=C CERUNWI*ETOLEN*£TOWID>METOLEN**2.)* 

2  SQRTtET0LEN**2.*ET0WID**2.-ERUNWI**2.)>/ 

Z  (ET0LEN**2.«ET0WID«*2.) 

CRUNL E=ER UNLE - < 6. • STAGE V ) 

BETL£NsSRL-t6.*STA0EV) 

CUTSs  tCRUNLE/BETLEN >-l • 

NCUTSsINT (CUTS  > 

REMsCUTS-NCUTS 

.  IFtRER  «GT •  0.011  NCUTSsNCUTS*l 
PCSTARsPCLOSE •«  <1 ./NCUT$> 

SPACE s t ERUNLE-ETOLEN) /(NCUTS) 

0 VRLAPs(ETOLE N -SPACE) /2. 


V 


CUTLOCI 1 ) =OVRL AP+SPACE-0  AHR AD 
00  100  I=2tNCUTS 

CUTLOCI I)=CUTLOC (I -1 )*SPACE 
100  CONTINUE 
RETURN 
ELD 

SUBROUTINE  BOUNOS(MIN«MAX) 

OIHENSION  NUMB0MI2Q )« A1HLOC(20) 

COMMON  /INPUT/  RUNL£N«RUNUlD«TOLEN«TOUIOtUtCEP 

COMMON  /TRANS/  ERUNLE«ERUNMI«ETOL£N»ETOUIOtDAHR  AO i ST ADC V 

COMMON  /GRP/  NSPWt STEP 

COMMON  /RUAT/  SRL 'PCSTAR .NCUTS.CUTL0CC20 > 

COMMON  SIMCAL 
LOGICAL  SIMCAL 
C 

MIN-ERUNUI/ETOUID 

C 

NUMAIHsRIN 

SPACE sRUNUIO/ IMIN*1 ) 

DO  20  11=1'20 
MAX=0 

00  10  I2sl«NUMAIH 
NUM80H4 121=11 
AIHLOCI 12 )=I2*SPACE 
MAX=MAX*I1 
10  CONTINUE 

IF ( NUHAIM  «E0.  1  .AND.  NUNROH(l)  .EG.  11  THEN 
CALL  PR  ONE  (  NUMA  I M  #  NUMBOM »  AI MLOC • PC) 

IF  (PC  .GE.  PCSTAR)  GO  TO  90 
ELSE  IF( SIMCAL)  THEN 

CALL  SIM(NUNAIM*NUHBOM*AlHLOC«PC) 

IF(PC>GE* PCSTAR ♦•01 )  GO  TO  90 
ELSE  IF («NOT •SIMCAL)  THEN 

CALL  UN I ONI NUMA I M'NUHBOM • A I MLOC «PC ) 

IF  (PC  «6E.  PCSTAR)  60  TO  90 
ENO  IF 

20  CONTINUE 

90  RETURN 

CNO 

SUBROUTINE  SEARCH  (MIN »HAX  *NUHB »PC ) 

DIMENSION  AIM(20)«AINO(20)fAIMI f 20 ) t AIHLOCI 20 ) (NUMBOM ( 20 ) 
DIMENSION  LAIM(S)«IAIM(S) 

OIHENSION  N8MINT(20)*AIM1NT(20) 

COMMON  /INPUT /  AUNLEN.RUNBID»TOLEN.TOMIO.M»CEP 

CCMMON  /TRANS/  ERUNLC 'ERUNU1 »ETOLEN«ETOWIOtOAMR AO* STAOCV 

COMMON  /6pP/  NSPNtSTCP 

COMMON  /RUAY/  SRL #PCS TAR.NCUTS »CUTLOC (20 ) 

COMMON  /BEST/  FCBCSTI20 ) »NALBST(20 ) .NBMBST (20 »5) . AIMBSTI 20 *5 ) 
COMMON  SIMCAL 
LOGICAL  SIMCAL 

C 

CNTR3 RUNWl 0/2 • 


LCNTRsCNTlUl. 

T  CL=0 *0 

IFCSIHCAL)  TOL-0«01 
C 
C 

DO  900  LOOP  =  PIN.  MAX 
NUHBsLOOP 
C 

100  NUHAIH=1 

NUHB0HC1)=NUHB 
AIHLOCC 1 > =CNTR 
IF  (NUMB  .EQ»  1)  THEN 

CALL  PRONE < NUNAI M .NUHBOH* AI HLOC . PCBEST Cl») 
NALBSTC1 )=NUHA !M 
NBHBST  Cl .1 ) =NUHBCMC1 ) 

AI MBSTC1 .1 l-AIHLOCC 1 ) 

GO  TO  800 
END  IF 

IFCSIHCAL)  THEN 

CALL  SI  H  CNUHAI  H«  NUHBOH.  AInHLOC*  PCBEST (NUHB)  ) 
ELSE 

CALL  UNION C NUHA1H. NUHBOH. AIHLOC. PCBEST C NUHB ) ) 
END  IF 

NALBSTCNUHBIsNUHAlH 
NBHBST  C NUHB .1 )=NUHB0HC1 ) 

AIHBST CNUHB.l )sAIHLOCCl ) 

C 

200  NUHAIHS2 

NREH=  MOOC  NUHB . 2 ) 

IFCNREH  *EQ.  1 )  GO  TO  300 

PC2=0 .0 

AIHCDsO.O 

AIHC2  IsRUf.UIO 

NUHB0HC1 )=NUHB/2« 

NUHB0HC2 ) -NUHB/2* 

C 

00  230  Jl'l.LCNTR.S 
AI HL0CC1 )= Jl-1 • 

A I HL0CC2 )-RUNUIQ-C  Jl-1 • 1 
IFCSIHCAL)  THEN 

CALL  SI H CNUHAIH . NUHBOH* AI HLOC .PCN ) 

ELSE 

CALL  UN IONCNUHAIH. NUHBOH. AIHLOC.PCN) 

END  IF 

IFCPCN  .GT.  PC2)  GO  TO  220 
LAST  sJl-10 
PC2-0 .0 

DO  210  J2sLAST.Jl»l 
AIHL0CC1 )=J2 
AIHLOCC  2 ) -RUNUI D-J2 
IFCSIHCAL)  THEN 

CALL  SIHCNUHAIH. NUHBOH. AIHLOC. PCN) 

ELSE 

CALL  UNICNCNUH A IH.NUHBOH. AIHLOC. PCN) 
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END  IF 

.  IFCPCN  .LE.  PC2 )  60  TO  240 
PC2=PCN 

AIM(1)=AIMLCC(1) 

AIM(2)=AIML0C<2> 

210  CONTINUE 

220  PC2=PCN 

AIMCl)sAINLOC(l) 

AIM(2)=AIHL0C(2> 

230  CONTINUE 

240  0IF=PC2-PCBEST(NUHB) 

IFfDIF.GT.TOL)  THEN 
PC BEST (NUMB) SPC2 
NALBST (NUMB)=NUMAIH 
NBMBST (NUHBvl ) =NUHBOM(l ) 

NBHBST (NUMB* 2)=NUMB0M(2 ) 

AI MBST (NUMB* 1 ) =AIM(1 ) 

A I MBST (NUMB *2 )=AIM(2) 

ELSE 

GO  TO  800 
END  IF 

IFCNUHB  .EO.  2)  60  TO  800 
C 

300  NUMAIHS3 

LBOHB=(NUMB-l >/2 
PC3=0 .0 
PCINT  =0 .0 

DO  360  NB0H=1 *  LBOHB 
NUHBOM (1 )=NB0M 
NUMBOM (2 ) -NUMB-2  *NBON 
NUHB0H(3)=NB0M 
AIM (1  >  =  0.0 
AIM(2 )=CNTR 
AIM<3 )=RUNU1D 
00  330  J1=1*LCNTR*5 
AIFL0CC1 >=U1-1. 

AIHL0C(2)=CNTR 

AI BLOC ( 3)=RUNMID-( Jl-1 I 

IF (SIMCAL)  THEN 

CALL  SIHtNUMAIM* NUHBOM* AlMLOCtPCN) 
ELSE 

CALL  UN ION(NUMAIM*NUHBOH»AIHLOC*PCN) 
END  IF 
OIF=PCN-PC3 

IF(DIF.GT.TOL)  60  TO  320 
LASTsJl-10. 

PC  3-0«0 

00  310  J2=LAST«J1*1 
AlMLOCd  )=J2 
AIMLOCt 2 ) =CNTR 
AIHLOCC  3 )=RUNWID- J2 
IF(SIMCAL)  THEN 

CALL  S1M(NUHAXM*NUMB0M»AIML0C*PCN> 
ELSE 


CALL  UN I ON < MJMA IN* NUHBOM  * A I ML OC *PCN I 
.  ENO  IF 

OIFsPCN-f C3 

IF  (DIF.LE.TOL)  60  TO  3*0 
PC3=PCN 

AIM(l)=AIMLCCtl> 

AIH<2)=AIHL0C(2) 

AIM(3)=AIMLOC(3> 

CONTINUE 

PC3=PCN 

A I M(l)=AIMLOC(l ) 

AIM(2)=AIML0C(2> 

AIH(3)-AIHL0C(3) 

CONTI NUE 

IF<PC3.GT.PCINT>  THEN 
PCINT=PC3 
NALINT=NUMAIM 
DO  350  1=1 *3 

NBMINT I I >=NUMBOMCI > 

AI MINT (I)-AIM(I) 

CONTINUE 
END  IF 
CONTINUE 

I FC  PC  I NT .6T. PC BEST (NUMB) )  THEN 
PC REST (NUMB)=PCINT 
NALBST (NUMB )=N ALINT 
DO  370  1=1 • 3 

NBMBSTt  NUMB* I )=NBMINT (I ) 
AIMBST(NUMB»I)=AIMINT(I) 

CONTINUE 

ELSE 

GO  TO  800 
END  IF 

IFCNUMB  .EQ.  3)  GO  TO  800 

NUMAI M=4 
NKEM= MOD (NUMB *4 ) 

IFCNREM  .EQ.  1  .OR.  NREM  .EQ.  31  GO  TO  500 

PC4=0.0 

PCINTsO.O 

LB0M8=(NUMB-2>/2 

DO  493  NBM=1 tLBOMB 

SPACE =R UN MI 0/ (NUMAIM+1 ) 

NUMBOMC 1 ) =NBM 
NUMBOMC  4) =NBM 
LBMB2  =  < NUHB-2  *  NBM I /2 
DO  493  NBM2=1 *LBMB2 

NUMBCM<2  >=NUMB-2*NBM 
NUHBOM (3  >  =  NUMB-2*NBM 
DC  410  I=1»NUHAIM 
AIMI)  =  I*SPACE 
LAIMd  )-AIM(I) 

CONTINUE 
IF(SIMCAL)  THEN 


UUU*  u  u  u 


CALL  SININUMAIM»NUHBOM* AIM»PC4 ) 

ELSE 

CALL  UNION4NUMAIH*WUNBON.AIM#PC4> 

END  IF 

♦30  STEP1=5. 

AIMI( 1)=AIHC1 ) 

AIMIC 4)=AIM 14 > 

AIHO( 1 )=AIH(1 ) 

AIHOC4l=AIM<4l 

HOVE  INNER  PAIR  OF  AIN  POINTS  IN  OR  OUT 

♦0  AIHI(2)=AIM(2)+STEP1 
AIHI(3):AIN(3l-SrEPl 
IF4SIMCAL)  THEN 

CALL  SIM<NUMA1M«NUHB0M»AIMX »PCI ) 

ELSE 

CALL  UNI0N(NUMAIR«NUM80M(AIMI«PCI ) 

E»D  IF 

AIH0l2|sAIN(2)-SrEPl 
AIHO(3)=AIM(3) ♦STEP1 
IFCSIMCALl  THEN 

CALL  SIN<NUMAIM«NUMBOM*AIMO«PCO> 

ELSE 

CALL  UNION (NUHA1 H* NUMBOH» AIHO.PCO > 

ENO  IF 
0IF1=PC0-PC4 
DIF2=PC4-PCI 
DIF3=PCI-PC4 
DIF4=PC4-PC0 

IF4DIF1.GT.TOL  *AND.  DIF2.GT.T0L>  THEN 
PC4=PC0 
AI N(2)sAIH0(2> 

AINC3>sAIH0(3» 

ELSE  IFC0IF3.GT.T0L  .AND.  DIF4.GT.T0L>  THEN 
FC4sPCI 

AIM<2>=AIMI <2) 

AINC3I2AIHIC3) 

ELSE  IF(0IF4.GT.T0L  .AND.  0IF2.GT.T0L>  THEN 
IFISTEP1  -LE.  0.5>  GO  TO  450 
STEPlsSTEPl/2. 

GO  TO  440 
E»D  IF 
• 

MOVE  OUTER  PAIR  OF  AIN  POINTS  I,N  OR  OUT 

450  AIMI(2>SAIM<2> 

AINI(3>*AIM<3> 

AIM0<2)sA1H<2> 

AIH0(3>sAIN(3> 

S*  EP2=5. 

460  AIMI<l>sAXHU)*ST£P2 
AIMI(4»3AIN(4)*STEP2 
IFtSIMCAL)  THEN 


CALL  SIHCNUMAIHtNUMBOMtAIHI ,PCI ) 

ELSE 

CALL  UMONC  NUHAIM.NUMBONt  AI MX  *PC I ) 

ELO  IF 

AIMOCl>=AIM(l I-STEP2 
AIM0C4>=A1H(4>»STEP2 
IFCSIPCAL)  THEN 

CALL  SIMCNUMAIH»NUHBOH»AIHO.PCO> 

ELSE 

CALL  UNION < NUHAI MtNUHBOHt AIHOtPCO ) 

END  IF 

DIF1=PC0-PC4 

OIF2=PC4-PCI 

DIF3=PCt-PC4 

DIF4=PC4-PC0 

IFC0IFl.6T.r0L  .AND.  0XF2.GT.T0L>  THEN 
PC4=PC0 
AIMC1>=AIMCC1> 

AIH(4>=AIM0C4) 

ELSE  IFC0IF3.GT.T0L  .AND.  DIF4.GT.T0L>  THEN 
PC4=PCI 

AINCl>sAIHIIl> 

AINC4>  =AIHI C4> 

ELSE  IFC0IF4.GT.T0L  .AND.  OIF2.GT.TOL>  THEN 
IF  CSTEP2 .LE.0.51  GO  TO  470 
STEP2sSTEP2/2« 

GO  TO  460 
END  IF 
C 

C  ROUND  AIN  POINT  LOCATIONS  TO  NEAREST  FOOT. 

C 

4  70  CONTINUE 

DO  480  1  =  1*4 

I'  .NCI >-AIM< I > ♦O .5 
480  CONTINUE 
C 

C  CHECK  MOVEMENT  FROM  LAST  ADJUSTMENT. 

C 

IFCLAIMC1 >.E0.IAIM(1)  .AND.  LAIM(2).EQ.IAIMC2> >  GO  TO  490 
00  485  I =1*4 

AXMCI>-I AIM(I> 

LAIMCDsIAIMCI) 

485  CONTINUE 
60  TO  430 

490  IFCPC4.GT.PCINT)  THEN 
PCINT=PC4 
NALINTsNUMAIM 
DO  492  1=1*4 

NBHINTf  I  >=NUMBOMCI  >’ 

AIRINTC I >=AI MCI ) 

492  CONTINUE 
ENO  IF 

493  CONTINUE 

IFCPC INT.GT .PC BEST (NUMB>  >  THEN 
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non 


PC0EST(NUM8>=PCINT 
NALBST(NUMB  >=NALINT 
'  00  495  1=1*4 

, NBM8ST ( NUMB* I )=NBMINT(I) 
AIMBST(NUMB«I)=AIMINT(I> 

495  CONTINUE 

ELSE 

GO  TO  BOO 
END  IF 

IF(NUMB.EQ.4)  60  TO  000 
C 

500  NUHAI M=5 

PC5=0.0 
PCI  NT  =0.0 
LBOM8=(NUMB-3>/2 
00  593  NBM=1*LBQMB 
SPACE =RUNUI0/CNUMAIH»1) 

NUHB0MC1 >=NBH 
NUMBOMt 5) =NBM 
L8M2=(NLMB-2*NBH)/2 
00  593  NBM2=1«LBM2 
NUMB0H(2)=NBM2 
NUMBOM(3)=NUMB-2*NBM-2*NBM2 
NUMB0MC4 )=NBM2 
DO  510  1=1* NU MAIM 
AIMUI=I*SPACE 
LAIMTI )=AI M(I ) 

510  CONTINUE 

IF  (SIMCAD  THEN 

CALL  SIMTNUMAI M«NUMB0M* AIM*PC5 ) 
ELSE 

CALL  UNIONTNUMAI M*NUMB0H* AIM «PC5 ) 
END  IF 

530  S~EP1=5. 

aimki)=a:m(1) 

AIMIC 3)=CNTR 
AIMI ( 5) =AIM (5 ) 

AIM0C1  )=AIM  (1  ) 

AIM0( 3)=CNTR 
A1MCT  5)=AIM(5) 

MOVE  INNER  PAIR  OF  AIM  POINTS  IN  OR  OUT. 

540  AIMH  2)  =  AIM(2)>STEP1 

AIHI<4)=AIMC4>— STEP1 
IF  (SIMCAD  THEN 

CALL  SIM(NUMAIM»NUHBOH»AlMItPCI> 
ELSE 

CALL  UNION (NUNAIM.NUMBON* AI MI *PCI) 
END  IF 

AIM0(2)=AIM(21-STEP1 
AIM0(4)=AIM(4) ♦STEP1 
IF  (SIMCAD  THEN 

CALL  SIM(NUMAIM*NUMBOM*AIMO.PCO> 


98 


non 


ELSE 

CALL  UNlONtNUMAIM»NUMBOM*AIHO#PCO> 

E>D  IF 
0IFlsPC0-PC5 
DIF2=PC5-PCI 
OIF3=PCI-PC5 
OIF4=PCE-PCO 

IFCOIF1.GT.TOL  .AND.  DIF2. GT.TOL)  THEN 
PC5=PC0 
AI M (2) =AIM0 <2 > 

AIM(4>=AIM0(4> 

ELSE  IFfOIF3. GT.TOL  .ANO.  0IF4.GT .TOL>  THEN 
PC5sPCI 

AIHC2MAIMII2) 

AINC4) =AINI f4> 

ELSE  IFf0IF4.6T.T0L  .AND.  DIF2.6T.T0L>  THEN 
IF  tSTEPL.LE.Q .5)  60  TO  550 
STEPl=STEPl/2. 

60  TO  540 
END  IF 

MOVE  OUTER  PAIR  OF  AIM  POINTS  IN  OR  OUT. 

550  AIMI<2>sA!M<2> 

A IM I ( 4 >  =  AIM  14 ) 

AIM0(2)-A1M(2> 

AIMOf  4>sAIM(4 ) 

STZP2s5. 

560  AIMIfllsAIMfl > ♦STEP2 
AIMIf 5>sAIMC5>-STEP2 
IF  fSIMCAL)  THEN 

CALL  SIM(NUHAIM«NUHBOM*AIMI tPCI > 

ELSE 

CALL  UNION! NUMAIMtNUMBOMf AIMI *PCI > 

END  IF 

AIMOf l)sAIMfl>-STEP2 
AIMOf 5>=AlHf5)*STEP2 
IF  fSIMCAL)  THEN 

CALL  SIMfNUMAI M»NUMBOM»AIMOtPCO ) 

ELSE 

CALL  UNIONf NUMAIMf NUMBOMt AI HO »PCO ) 

END  IF 

DIFlsPCO-FCS 

DIF2sPC5-PCI 

DIF3=PCI-PC5 

OIF4SPC5-PCO 

IFfDI FI .GT.TOL  .ANO.  DIF2.6T.T0L)  THEN 
PC5*PC0 

AlMfl)=AIMOfl) 

AIHfS)=AIH0(5> 

ELSE  IFfDIF3. GT.TOL  .AND.  DIF4.6T.T0L>  THEN 
PC5*PCI 

AI  Mf  1  >=AIM  HI  > 

AIMf5)-AIMIfS) 
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ELSE  IFCDIF4.GT.T0L  .AND.  OIF2.Gf.TOU  THEN 
IF (STEP2.LE.0.5)  GO  TO  570 
STEP2=STEP2/2. 

GO  TO  560 
ENO  IF 
C 

C  POUND  OFF  AIM  POINT  LOCATIONS  TO  NEAREST  FOOT. 

C 

570  OC  580  1  =  1*5 

IAIHC I)=AIM(t >*0.5 
580  CONTINUE 
C 

C  CHECK  MOVEMENT  FROM  LAST  ADJUSTMENT. 

C 

XFCLAIMfl).EQ.tAIM(l)  .AND.  LAIN! 2 I.CQ.I AXM(2> )  GO  TO  590 
OC  585  1=1*5 

AIMCI )=I AIM(I ) 

LAIMCI )=IAIMII> 

585  CONTINUE 
GO  TO  530 

590  XF(PC5.GT.PCINT)  THEN 
PC  I NT=PC5 
NALINT=NUMAIM 
DQ  592  1*1*5 

NBMINTl 1 )=NUMBOMC l) 

AIMINT(I)=AIM(I > 

592  CONTINUE 
ENO  IF 

593  CONTINUE 

IFCPC INT.GT .PCBEST (NUMB ) >  THEN 
PC BEST (NUMB )=PCINT 
NALBST (NUMB)=NALINT 
00  595  1=1*5 

NBMBSK  NUMB*  I  )*NBMINT (I > 

AIHBST ( NUMB* I )=AIMINT <11 
595  CONTINUE 

ENO  IF 
C 
C 

800  0IF=PCBEST<NUMB1-PCST  AR 

IF(OIF.GE.TOL)  GO  TO  1000 
C 

900  CONTINUE 
C 

C  FOUNO  THE  8EST  STRATEGY  TO  CLOSE  THE  RUNUAY. 

C 

1000  PC*PCB£ST (NUMB) 

RETURN 

E»,D 

SUBROUTINE  RESETS (NPERCT * NUMBER *PCLEST ) 

DIMENSION  NPERCT <20 ) 

COMMON  /INPUT/  RUNLEN«RUNUIO»TOLEN*TOUlD*U*CEP 
COMMON  /KUAY/  SRL .PCSTAR *NCUTS*CUTLOC C20 » 


COMMON  /BEST/  PCBEST120  )#NALBST  (20 ItNBMBST  !20  #5) •  AIHBST!20« 5) 
C 

PRINT  10  «RUNLENtRUNWID 
PRINT  2OtTOLEN»TOWI0 
PRINT  30 
PRINT  40(W 
PRINT  SO* CEP 
PRINT  60 tPCLEST 
PRINT  70 • NUMBER 
PRINT  80 
PRINT  90 

DO  200  IsltNCUTS 
NUMBs NPERCT d 1 

00  200  d=l tNALBST ( NPERCT (1)1 

PRINT  100t CUTLOC d  ) « AIMBST !NUMBtJ) • NB MOST ! NUMB «d) 

200  CONTI NUE 

10  FORMAT! *1 RUNWAY* tl3XtF8«0t*  BY  **F8.0f*  FEET* I 

20  FORMAT(*OMIN  LAUNCH  UINOOU  **F8.0t*  BY  *«FB.Ot*  FEET • 1 

30  FORMA T  !  *0 WEAPON  CHARACTERISTICS* ) 

40  FORMAT! *0*t5X»* YIELD  *»F7#2**  POUNDS*) 

90  FORMA T ! *0 **5X # *CEP  '  **3X»F4«0t*  FEET*) 

60  FORMAT!*OPR OB ABILITY  OF  CLOSURE  *»F4.2) 

70  FORMAT! *OTOTAL  NUMBER  OF  WEAPONS  **I3) 

80  FORMAT! *0 AIM  POlNTS*«7Xt *NUMBER  OF*) 

90  FORMAT! *0!LEN6TH«  WIDTH) • «2X« *WEAPOMS* %t ) 

100  FORMAT!*  *tF7.2tF8.2«4XtI2) 

RETURN 

ENO 

SUBROUT INC  NUMSET INPERCT • NUMBER  tPCLEST ) 

DIMENSION  NPERCT  (.0) 

COMMON  /RWAV/  SRL tPCSTAR tNCUTSt CUTLOC C 20 ) 

COMMON  /BEST/  PCB£ST!20)vNALBST!20)tNBHBST!20«S)«A!MBST!20»S) 
C 

PCSTARsl.O 

NREM-MOO! NUMBER tNCUTS) 

NUMBs  NUMBER/NCUTS 
NUMBl*MiRB*l 
DO  100  IsltNCUTS 
NPERCT !I >sNUMB 
100  CONTINUE 

CALL  SEARCH!NUMB«NUMB1 tNUMBl »PC) 

DO  200  IslvNREM 

’ NPERCT II >s NPERCT (I )*1 
200  CONTINUE 

PCLESTsl.O 
00  300  IsltNCUTS 

PCLEST*PCLEST*PCBEST!NPERCT!I) ) 

300  CONTINUE 
RETURN 
ENO 

SUBROUTINE  INOEP 

01  HEN  SION  NUNBOM ! 20 ) t A I MLOC ! 20 ) t CUTLCI ! 20 ) 

COMMON  /RWAV/  SRLtPCSTAR tNCUTS .CUTLOC !20 ) 


COMMON  SIMCAL 
LOGICAL  SIMCAL 


PCtNO=1.0 
NUMB- 0 
PRINT  100 
PRINT  200 
PRINT  300 
PRINT  400 
READ* »I CUTS 

IF(ICUTS.LT.NCUTS)THEN 
PRINT  &00 .NCUTS 
PRINT  310.ICUTS 
RETURN 

ELSE  IF(ICUTS.GT.NCUTS)THEN 
PRINT  800.NCUTS 
PRINT  820 i ICUTS 
RETURN 
END  IF 

DO  20  1=1. ICUTS 

RE AO*«CUTLCI (I).NUMAIM 
PRINT  500 .CUTLCI (I ) 

IF  (l.EQ.ll  THEN 
OISTsCUTLCI < I ) 

ELSE 

OIST=CUT LCI C I )-CUTLCI (1-1 > 

END  IF 

XFCOIST.GT.SRL)  THEN 
PRINT  900.SRL 
RETURN 
END  IF 

DO  10  Jsl.MJNAIM  x 

PEAD*«AIMLOC!J) .NUMBOM(J) 

PRINT  60  O.A1MLOC! J)«NUMBOH!J) 

NUMB-NUMB'* NUMBOMC  J) 

CONTINUE 

IF(NUMR.EQ.1)THEN 

CALL  PR0NECNUMAIM.NUMB0M.AIML0C.PC) 

ELSE  IF  (SIMCAL)  THEN 

CALL  SIMCNUMAIM.NUMBOM.AIMLOC.PC) 

ELSE  IF  (.NOT.  SIMCAL)  THEN 

CALL  UNIONCNUMAI M.NUMBOH.AIMLOC.PC) 

END  IF 

PC I ND=PCINO*PC 
CCNTINUE 
PRINT  700 .PCI NO 

FORMAT ( #1 INDEPENDENT  ANALYSIS*) 

FORMAT <*0 INDEPENDENT  STRATEGY:*) 

FORMAT! *0CUT*  «T12.*AIM  POINT* .T24 . 'NUMBER  PER*) 
FORMAT! •  LOCATION  *  *T1 2 .* LOCATION*  *T24.*AIM  POINT*) 
FORMA r (*0*.F8«0) 

FORMAT! •  *.T14,F8«2.T27.I3) 

FOR MATC'O PR OB ABILITY  OF  RUNWAY  CLOSURE  IS  *.F4.2) 
FORMAT ( *0  •  ERROR  NUMBER  OF  CUTS  MUST  EQUAL  **I4) 


810  FORMAT!'  *,I4f*  CUTS  ARE  TOO  FEU  TO  6UARENTEE  RUNUAY  CLOSURE' > 

820  FORMAT!'  'iMi*  CUTS  ARE  MORE  THAN  THE  REQUIRED  MINIMUM  NUMBER  OF 

1CUTS') 

900  FORMAT  CO  ERROR  DISTANCE  BETWEEN  CUT  LOCATIONS  IS  TOO  LARGE  TO 

I  CLOSE  RUNUAY. •*/*•  THE  MAXIMUM  DISTANCE  IS'«F10.2»'  FEET  ALONG  TH 
IE  RUNUAY') 

RETURN 

END 

SUBROUTINE  EV ALINPERCT. NUMBER. PCLEST I 
DIMENSION  NPERCT (20  I 

CCMMON  /BEST/  PCHEST! 20 ) tNALBST (2  0 >  *NBNBST! 20  *5) « AIMBST 120 tS) 
COMMON  /RUAY/  SRL « PCS TARf NCUTSt CUTLOC !20 > 

C 

NUMBERsO 

PCLESTsl.O 

DO  100  I-ltNCUTS 

NUMBER=NUMBER»NPERCT( I) 

PCLEST sPCLEST*PCBEST!NPERCTlI>) 

100  CONTINUE 


RETURN 

ENO 
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